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CHAPTER  I 


BASIC  CONCEPTS 

1.1  Introduction,  Objective 

In  the  analysis  of  linear  dynamic  systems  the  mathematical  des¬ 
cription  can  be  formulated  as  a  single  ordinary  linear  differential 
equation  with  constant  coefficients.  Upon  application  of  transform 
methods  a  polynomial  commonly  called  the  "characteristic  equation"  is 
obtained.  An  essential  part  of  most  analysis  problems  is  the  evalua¬ 
tion  of  the  roots  (zeros)  of  this  characteristic  equation.  Mitrovic  (1,2) 
has  developed  a  graphical  technique  which  permits  ready  evaluation  of  the 
roots  as  a  function  of  the  values  of  two  of  the  coefficients  of  the  poly¬ 
nomial. 

While  there  are  many  methods  for  evaluating  the  roots  of  a  poly¬ 
nomial,  few  of  these  methods  are  of  significant  advantage  in  engineering 
synthesis  problems.  Mitrovic' s  method,  because  it  relates  the  root  values 
to  coefficient  values  in  a  convenient  graphical  representation,  provides 
a  very  useful  approach  to  design. 

Most  of  the  methods  available  for  analysis  and  design  are  either 
one  parameter  methods  (such  as  the  root  locus  method)  or  evaluate  dynamic 
characteristics  indirectly,  or  both  (such  as  most  frequency  domain  methods). 
Mitrovic 's  method  is  basically  a  two  parameter  method,  but  is  readily  ex¬ 
tended  to  three  parameters .  It  is  possible  to  extend  the  method  to  handle 
more  than  three  parameters,  but  the  practical  value  of  such  techniques 
has  not  yet  been  established. 

The  original  development  by  Mitrovic  (1,2)  considered  only  two  co¬ 
efficients  of  the  polynomial,  i.e.,  for  an  expression  of  the  form 
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n-1 


(1-1) 


Che  original  development  considers  analysis  and  synthesis  only  in  terms 
of  aQ  and  a^.  More  recently  (3,4)  the  basic  concepts  have  been  gener¬ 
alized  to  permit  analysis  and  design  in  terms  of  any  two  coefficients. 

The  objective  of  this  research  paper  is  to  provide  detailed  analyses, 
procedures,  equations  and  techniques  which  will  reduce  the  more  recent 
advances  to  a  more  practical  level.  This  is  to  be  done  by  providing: 

a)  detailed  stability  analyses  and  rules  for  applying  them  to  any 
of  the  Mitrovic  planes. 

b)  tabulated  equation  pairs  to  assist  the  computations. 

c)  procedures  and  techniques  for  sketching  the  curves  with  a  min¬ 
imum  amount  of  computational  labor. 

d)  tabulated  relationships  to  aid  the  sketching. 

e)  a  digital  computer  program  to  assist  with  detailed  computa¬ 


tions 


f)  some  selected  illustrations  of  the  results  in  design  applications. 
1-2.  Review  of  Mitrovic's  Method. 

If  all  of  the  roots  of  a  polynomial  be  inside  some  area  in  the 
s-plane,  then  proof  of  this  can  be  established  by  enclosing  the  area  by 
a  contour ,.  {napping  the  contour  onto  a  polar  plane  through  the  charact¬ 
eristic  polynomial  as  a  mapping  function,  and  analyzing  the  polar  contour 
with  the  Principle  of  Argument.  This  is  essentially  Mitrovic's  Method. 

He  chose  as  mapping  contours  only  the  imaginary  axis  or  a  radial  straight 
line  in  the  left  half  plane,  closing  his  contour  through  a  circular  arc 
of  infinite  radius  so  as  to  enclose  all  or  part  of  the  left  half  plane. 
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With  mapping  contour  defined  as  above,  a  point  on  the  contour  is  de¬ 


fined  by 


=  -w  sin  0  +  j  a)  cos  0 


n  n 


=  -  C  con+  j  u>n/  1  -  c2 


(1-2) 


Inserting  equation  1-2  in  1-1  and  requiring  that  reals  and  imaginaries  are 


zero  independently  provides  two  equations  in  terms  of  u)n,  the  coefficients, 


and  certain  functions  of  £.  But  £  is  a  known  number  and  for  the  usual  case 
of  analysis  the  coefficients  are  known.  Thus  two  equations  in  one  unknown 
(Wn)  are  obtained.  Hitrovic's  procedure  was  to  define  (select)  two  coeffi¬ 
cients  as  unknowns,  thus  providing  two  equations  in  two  unknowns  and  a 


common  parameter,  co  .  He  solved  the  equations  for  the  unknown  coefficients 
in  terms  of  uj^.  Then' choosing  a  value  for  £  and  letting  u)n  take  selected 


values,  curves  can  be  obtained  by  plotting  one  coefficient  as  ordinate  and 


the  second  as  abscissa  with  u)^  as  a  parameter. 


If  this  is  done  for  a  polynomial  for  which  all  coefficients  are 
known  numerically,  then  the  actual  values  of  the  coefficients  which  were 
designated  variables  are  the  coordinates  of  a  single  point  on  the  Mitrovic 
plot.  Mitrovic  showed  how  to  evaluate  stability  and  all  left  half  plane 
roots  of  the  polynomial  from  his  curves  and  the  location  of  this  one  point. 
If  this  point,  called  the  M-point,  is  moved  to  a  new  location  new  roots  and 
new  coefficient  values  are  defined.  Then  the  physical  system  can  be  provi¬ 
ded  with  the  dynamics  specified  by  the  new  roots  if  the  system  can  be 
changed  so  that  the  two  designated  coefficients  assume  their  new  values 
without  changing  any  other  coefficients.  This,  in  essence,  is  the  result 
provided  by  Mitrovic,  and  to  it  he  added  considerable  detail  regarding 
specific  techniques  for  using  the  method  to  design  feedback  control  systems. 
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In  his  treatment  Mitrovic  defined  as  variables  only  the  coefficients 


a^  and  aQ  (see  equation  1-1).  Elliott,  Heseltine  and  Thaler  (3)  used  sev¬ 
eral  other  coefficients  as  variables  in  specific  problems  involving  the  de¬ 
sign  of  feedback  compensation  for  control  systems,  but  did  not  generalize 
these  results.  Later  Siljak  made  an  independent  study  which  generalized 
Mitrovic' s  results  to  any  pair  of  coefficients,  but  did  not  expand  or  apply 
his  results.  This,  then,  was  essentially  the  state  of  development  of 
Mitrovic's  method  when  this  research  was  undertaken;  the  basic  method  using 
only  the  aQ  and  a^  coefficients  is  well  developed  and  techniques  are  avail¬ 
able  for  applying  it  to  many  design  problems.  Mitrovic  (1)  and  Siljak  (5) 
also  applied  this  method  to  sampled  data  systems.  The  basic  theory  of  the 
generalized  method  has  been  established,  and  a  few  techniques  have  been 
developed  to  apply  it  to  feedback  compensation  design.  In  general,  practical 
ways  to  interpret  stability,  and  to  apply  the  generalized  theory  to  analysis 
and  design  are  lacking.  Some  advances  in  these  areas  are  presented  in  this 
report . 
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CHAPTER  II 


THE  GENERALIZED  METHOD 


2.1  Siljak's  Results. 

Through  a  formal  derivation  and  rigorous  proof  Siljak  obtained,  in 

symbolic  form,  a  generalized  expression  for  the  Mitrovic  Curves  in  terms 

of  any  two  selected  coefficients.  These  results  are;  for  Bi  =  Bj,  (i  >  j) 
"N  “*  “  _r 


Bi  = 


AA  “n  Vl  <C> 


B<J  = 


N 


E  Ak  ^n  ^k-i  (C) 

k  =0 


(C) 


where  A^  is  the  coefficient  of  the  k'th  order  term  of  the  original  polynomial, 
and <p^(C)  are  the  functions  defined  by  Mitrovic. 

2.2  The  Determinantal  Approach.  Tabulation  of  Equation  Pairs. 

Mitrovic' s  original  result,  (rearranged  as  in  Ref.  2)  expressed  with 


all  coefficients  designated  as  a's  is: 


a  =  -oo  [a, 
o  n  z 


+  a.  u L  CP,  +  a,  o  < p,+  ...  +  ar 


wn"2<Pn-l 


al=  a2<p2  con+  a3cp3  J+  a  4  +  *n<Pn  oj"1. 

where  the<p's  are  functions  of  £  and  are  tabulated  elsewhere  in  this  report. 

In  the  original  method  aQ  and  a^  were  designated  variables.  In  gen¬ 
eral  any  two  coefficients  may  be  chosen.  To  obtain  specific  equations  for 
any  two  coefficients  rearrange  the  above  equations  with  both  of  the  chosen 
coefficients  on  the  left  hand  side;  for  example  a2  and  ay 

a2  wn  *1  +  a3  “n  * 2  =  -ao  *  “n  [a4  “n  ^  +  —  +  a»  “"’2  Vl] 

2  3  n*  1 

-a0  <p  a)  -  a,  <p,  UL  =  -a.  +  a  <p.  ojI  +  .  .  .  +  a  <p  o 
2  2  n  3  3  n  1  4  4  n  n  n  n 
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These  equations  are  simply  simultaneous  equations  in  a?  and  a 

3* 

may  be  solved  by  the  usual  determinantal  methods,  i.e., 

a2  ^ 


a2 


A  *  a3  ” 


where  A  - 


6t)2  p- 

n  1 


'"n'Pz 


“n^2 

-o?  <f> 

n  3 


-a 


-oi  [a  to2<p.  +  ...+  a  ton"2<p  ,  1  UJ3  <p 
n  L  4  n  3  n  n  n-1  J  n  ^ 


-a.  +  a.  <p.  uP  + 
1  4  ^4  n 


•  +  a  <Pn  tiL 
n  n  n 


n-1 


-to2  <p_ 
n  3 


^3  ~ 


“n^l 


-to  <Pn 
n  2 


-a 


2  r  2  n-2  1 

-to  a.  to  <p_  +  . . .  a  to  <p  , 
n  L  4  n  3  nn  n-1  J 


-al  +  “4  ”4  “n  +  •••  + 


If  the  indicated  algebra  is  performed  the  desired  relationships  are  ob¬ 
tained,  as  tabulated  in  Table  2-1  which  follows; 


and 
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TABLE  I 


(L)  Bq  -  Bt 

Bd  ‘  [  A2  <4*1  «>  +  A3  <4*2  «>  +  •••  +  An  <4Vl  «>  ] 

B1  _  A2  “n  *2  «>  +  a3  <4  *3  K>*  ■■■*  An  “V  *n  «> 

(2)  B0  ■  Bz* 

Bo  *  fl1  «>  [+  A1  <S,*1  «>  +  A3  <4*1  «>  +  A4  <4*2  (°  + 

",An  %*n-2  ^  ] 

B2  *  % 2*2L  <i>  [  *A1  “n*l  «)  -  A3  <4 *3  «>  '  A4  <4*4  «>’ 

’  '  ’  "  Ati  “n  *n  ] 

4>4- 

(3)  B0  -  B3 

B0  -  «)  [-  un«>2  «)  -  a2  UB<P1  «>  +  A4  ^  (O  +  , 

...  +  An  wn'Pn_3  (C)  ] 

b3  ■  -“iiS1  «)  [  +  *!  «  9j«>  +  a2  <4<p2«>  +  A4  J<ek  CO  + 

•••  +  A„<4*„  «>  ] 

(4)  B„  -  B4* 

Bg  =  (Pi1  CD  [ax  Wn  cp3  (D  +  A2  o£  <P2'  (°  +  A3  wn  ^1  <C)  * 

A5<4*1  «>  -  *n“fiV4  «>  ^ 
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B4  "  ~un  [  +  A1  “n^l  K>  +  A2  "n^  (C)  +  A3  <4*3  (C)  + 

*5  <4*5  (O  +  -  +  A„  “On  «>  ] 


(5)  Bj  -  B2 


h  ’  '  %l  [  -  Ao«>2  «>  +  A3  <4*1  «>  +  a4  <4*2  «>  +  ••• 

+A„“>u-2  «>] 

B2  -  %  [‘  Ao*l  «>  +  a3  <4*2  «>  +  -+  A„  <4*n-l  <«  ] 


(6)  Bj  -  B3 


B1  *  “d1  vl1  «>  [-  A0*3  «>  •  A2  <4*1  «)  +  A4  <4*1  «>  + 

+  A»“n*n-3«>] 

«>[Ao+A2<4*1  «>+A4<4*3  «)+--+A»“>„-l 


<?>  Bj  -  B4 


Alt 


Bi  =  w;1  (P*1  «)  [-a0<p4  <t)  -  a2«>2( o  -  a3  <4^(0+  a5 

<4*1  «>  +  +  V4*n-4  «)  ] 

B4=-“n4  *3*  «>  LAo  *  a2  <4  *1  «>  +  a3  <4  *2  «>  +  A5  “4  *4«>  + 

'••  +  Ati  <4*n-l  K>  ] 
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(8)  B2  -  B3 


B2  =-u>;2  [-  A 003  (?)  -  Aj  Wn<P2  (O  +  A4  <4^  «)  +  A5  <4<P2  (O  ♦ 
•••+An"n%-3  «>  ] 

B3  *  “£3  [-  ko'l>2  «>  ‘  A1  wn  *^1  <°  +  A4  “n”2  «>  +  •••+  An 

<o  ] 

(9)  B2  •  B^* 

b2  -  -;2 ^2l  (O  [-  ao„4  «)  -  Aj  %<I>3  (O  -  a3  <4^  (O  + 

A5  %  ^1  +  •  •  •  +  An«>n.4  (O  ] 

»4~  “n4^  <»  [  -Ao^  «>  '  Al  <Si«,l  +  A3  «>  + 

*5  <4  ”3  <C+  -+An“SV2  «>  ] 

(10)  b3  -  B4 

b3  =  ^  v‘1  «>  [-  a04>4  (O  -  Aj  «n<p3  <0  -  Aju£^  (C)  +  a5  0)B  <dx 

«>+  •••+A„"X-4«>] 

b4  '“i4^'!  «>  [-  Ao»>3  «>  •  A1  %*>2  «>  -  a2  “n^l  «>  +  Vn 

*>2  «>  +  -  +An  “£  V3  «>  ] 

Note:  Those  with  *  and  **  have  singularities  at  £  =  0  and 
C  =  .5  respectively. 
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2.3  Discussion. 

The  equations  of  Table  I  are  all  of  essentially  the  same  form, 
and  may  be  plotted  or  sketched  in  essentially  the  same  way.  Inter¬ 
pretation  of  the  results  is  not  obvious,  however.  Certainly  if  the 
imaginary  axis  is  mapped  (£-  -  0)  it  should  be  possible  to  define  stab¬ 
ility  from  the  location  of  the  M-point,  but  stability  has  been  defined 
pictorially  only  for  the  Bq  vs  B.^  case.  Thus  the  interpretation  of 
stability  and  other  relaLed  topics  must  be  investigated,  and  this  is 
reported  in  Chapter  3. 

One  immediate  difficulty  appears  in  the  equations  for  even- 

even  or  odd-odd  coefficients,  such  as  B  -  B,;  B-  -  B,  ;  B,  -  B-5 

o  1  2  4  1  3 

Bq  -  B^.  All  such  equation  pairs  contain  a  factor  of  the  f orm <Peven  (C  )• 
Since  all  such  factors  contain  £  as  a  factor,  then  (peJen  (£  =0)  -*  <»  and 
the  equations  are  uninterpretable.  We  have  not  as  yet  found  a  fundamental 
explanation  of  this  nor  a  convenient  way  out  of  the  dilemma.  However, 
it  is  quite  certain  that  the  curve  defined  by  these  equations  is  a  finite 
curve,  for  example  any  selected  point  on  the  Bq  vs  Bn  curve  for  £  =  0 
is  also  one  finite  point  on  a  Bq  vs  B7  curve  for  which  has  the  de¬ 
signated  value.  Further  evidence  of  this  was  obtained  as  follows: 
equations  for  negative  C’  were  derived  and  Bq  vs  B2  Mitrovic  curves 
were  calculated  by  digital  computer  for  a  third  order  polynomial  using 
C  =  +  .003  and  +  .001.  The  results  are  shown  on  Fig.  2-1.  The  grouping 
of  the  curves  is  clear  evidence  that  the  B  vs  B„  curve  for  C  =  0  must 
lie  between  the  £  =  +  .001  curves,  and  thus  must  be  finite.  Tn  practice 
use  of  a  curve  for  (  =0.1  should  normally  be  satisfactory  for  engineering 
purposes . 

In  like  manner  for  certain  other  pairs  such  as  B  -  B-  and 
-1  0  3 

B^  -  B^  a  factor  of  (p^  (C)  obtained,  which  approaches  infinity  at 

C  =0.5  If  other  similar  singularities  exist  they  have  not  yet  been 

observed . 
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CHAPTER  XII 
STABILITY  CRITERIA 


3.1  Basic  Derivations 

A  stability  criterion  based  only  on  the  curve  of  on  the 

real  two  dimensional  plane  is  derived.  Any  interpretation  in  terms  of  a 
phase  plot  of  F  (s)  on  the  complex  plane  is  undesired,  because  the  former 
curve  is  plotted  first,  and  it  is  more  convenient  to  skip  the  phase  plotting, 
if  possible. 

Any  polynomial  F  (si  can  be  put  in  the  following  form: 


F  (si  =  f  (si  +  A, 


+  A. 
J 


(3-1) 


which  is  not  zero  if  s  is  not  a  root  of  F  (s).  Also  the  same  polynomial 


in  terms  of  Bi  and  is  given  by: 

F  (s)  =  f  (s)  +  Bi  s1  +  Bj  s'* 


(3-2) 

It  must  be  noticed  that  this  F  (s)  is  zero  for  any  value  of  s, 


provided  B..  and  satisfy  the  relations  given  already  in  the  Table. 

Thus  for  these  values  of  and  B ^ ,  equations  (3-1)  and  (3-2)  combine  to: 


F  (si  =  (A  -  B.)  s 
i  i 

Here  s  is  a  complex  number 

s  =  tu  /ff/  2  +  6 

Substitute  this  into  (3-3),  and 


+  (A.  -  B.)  s' 
J  J 


(3-3) 

(3-4) 


r(uln)  =(At  -  B.)  /l  (7T/  2  +  e>  +  (A.  -  Bj)  [j  (IT/  2  +  61 

=<*i  -  v  <4  V  <v  V  4  v  f3-5) 

e^  and  e^  are  unit  vectors  in  direction  of  — j  +  i  0  and  ^  +  j  0  re¬ 
spectively. 

Tt  is  evident  if  the  phase  plot  of  F(to  )  with  regard  to  to  varying 

n  si 

from  zero  to  infinity  encircles  the  origin  just  the  angle  of  n  V  /  2  +n  0  in 
the  counterclockwise  direction,  that  all  of  the  roots  have  a  damping  ratio  £ 
larger  than  £  =  sin  0 .  The  vector  relationship  in  equation  (3-5}can  be  used 
to  determine  some  basic  rules  for  the  interpretation  of  stability  on  the 
Mitrovic  plot*  These  rules  are  developed  in  the  following  paragraphs  for  most 
of  the  cases  listed  in  Table  I. 
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It  should  be  noted  that  stability  is  determined  by  mapping  the 
imaginary  axis  of  the  s -plane,  ie.  ,  by  the  Mitrovic  plot  for  £  —  0, 

A  rule  for  this  case  is  readily  developed.  When  designing  for  good  dynamic 
performance  the  location  of  the  m  -  point  can  also  determine  whether  a 
specification  £  =  is  satisfied.  Rules  for  this  are  also  developed. 

fl>  Bo  '  B1 

The  criterion  for  this  case  has  been  given  by  Mitrovic  for  both 
£  -  0  and  £  >  0  cases. 

It  must  be  noticed  that  this  belongs  to  the  special  group  of 

cases  where  the  rule  that  B  -  const,  must  be  cut  first  is  valid 

o 

even  for  £  >  0, 

Rules 

For  all  of  the  roots  to  have  a  damping  ratio  £  larger  than 
sin  § ,  the  Bq  -  const,  line  must  be  cut  first.  Cutting 
points  must  lie  alternately  on  JBc  =  const,  and  -  const, 
lines,  as  is  increased  from  zero  to  infinity. 

a)  bx  -  b2 

(i)  C  =  0 

This  leads,  0  being  zero,  to 
=  Ull’  e2  =  [n_, 
from  equation  (3-5). 

As  is  shown  in  Appendix  71,  it  is  valid  in  any  pair  that 

lim  F  (  ul  )  =  A  e  (3-6) 

Ui  n  c  o 

n  -*  o 

So  the  phase  plot  of  F  (  ui  )  starts  from  A  on  the  real  axis 

_  n  *  o 

at.  =  0,  and  must  cut  e,  first.  This  implies 

A2  =  B2  and  Al  >  Bj  (3-7) 

at  the  first  cutting,  as  can  be  seen  from  the  equation: 

F  (  a^  )  =  (  A2  -  B2  )  e2  +  (  Ar  -  )  uJn  ^.(3-8) 

This  is  illustrated  in  fig.  3-2;  the  provision  (3-7)  is  re¬ 
alized  when  B2  —  const,  line  is  cut  first. 


3-2 


(ii>  C  >  0 

In  this  case  each  unit  vector  becomes 
e1  =  /ff/2  +  6,  e2  =  fit  +  26. 

As  is  seen  from  fig.  3-3  of  the  phase  plot,  -e2  must  be  cut 
first.  This  implies  that 

A2<B2,  Ai=Bl  0-9) 

This  is  satisfied  only  when  =  const*  line  is  cut  first  as 
shown  in  the  same  figure* 

<3>  B2  -  B3 

(i)  c  =  o 

For  this  case,  from  equation  (3*5) t 
e2  =  ijL-  »  e3  =  /  3  tt/2  . 

Therefore  the  B2  =  const,  line  must  be  cut  first  as  is  ex¬ 
plained  in  fig.  3-4. 

(ii>  C  >  0 

Again  applying  equation  (3-5),  where 

e2  =/  77  +  2  6 ,  e3  =  )  3  77/  2  +  3  6  , 

it  is  easily  seen  the  point  cut  first  lies  on  the  direction 

m  — + 

vector  -  which  leads  to 

*  ^3  —  ®3  *  ^2  ^  ^2* 

Thus  it  is  concluded  that  namely  Bodd  must  be  cut  first  in 
this  case  also,  as  shown  in  fig*  3-5. 

<4>  E3  '  BA 

(i)  C  =  0 

By  the  same  analysis  as  the  previous  ones  it  is  evident  that 

the  B  s  const*  line  must  be  cut  first  in  order  to  have  a 
even 

stable  domain, 
fit)  C  >  0 

The  Bodd  s  const*  line  must  be  cut  first* 

(5)  B^  -  B- 
o  J 

For  this  case  equation  (3-5)  becomes 

F  (oo  )  —  (  A-  *  B_)  u?  e  +  (  A  -  B  )  e  . 
n  3  J  n  3  ooo 
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(i)  C  =  0 


-  /-3  lit.  2  ,  e  =  [o. 


(6)  Bj_  -  B4 


The  phase  plot  shows  -  must  be  cut  first,  which  Leads  to 

A,  <  B_  ,  A  =  B  . 

3  3  o  o 

This  implies,  as  is  .seen  from  fig.  3-6,  that  Bo  const,  line 
must  be  cut  first  in  B,  -  plane. 

(ii)  C  >  0 

®3  =  /3  77/  2  4-  3  6,  eQ  =  [o_. 


The  phase  plot  for  this  case  shows  that  the  point  cut  first  must 
satisfy 

A3<B3»  Ao5=Bo' 

which  indicates,  in  terms  of  BQ  -  Bj  curve,  that  the  Bq  * 
const,  line  must  be  cut  first  as  shown  in  fig.  3-7. 

Equation  (3-5)  becomes 

F(  wn  )  *=  (  A4  -  B4  )  e4  +  (Ax  wn  er 


(i)  C  =  0 


®4  =  Ll JL  =  /_0 _ ,  e1  =  j_nT2_ 


From  the  same  analysis,  it  is  evident  that  the  e^  line  must  be 
cut  first,  which  means,  in  terms  of  Bj  -  B,  plot,  B4  =  const, 
line  must  be  cut  first. 

(ii)  C  >  0  , 

®4  =  / 2  it  +  4  8  ,  =  /ff/  2  +  8 

which  indicates  that  the  =  const,  line  must  be  cut  first. 
All  of  the  remaining  cases  in  the  Table  have  singularities  at  £  =  0, 
but  this  difficulty  can  be  avoided  by  merely  starting  the  value  of  £  from  a 
very  small  quantity.  So  C  =  0  case  is  dropped  here  intentionally. 

(7)  B_  -  B~ 


£  is  always  assumed  to  be  larger  than  zero,  and  unit  vectors  in 
equation  (3-5)  becomes 

s  ts  ±  2  9>  ®n  =  /_Q_  • 
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Fig.  3-8  shows  that  -e2  must  be  cut  first,  which  means,  in 
Bq  -  B2  terminology,  Bq  -  const,  line  must  be  cut  first  in 
order  to  make  all  of  the  roots  have  £  larger  than  sin  0 . 

(8)  B  -  B, 

In  this  case  unit  vectors  of  equation  (3-5)  are 


Almost  the  same  analysis  as  before  leads  to  the  same  con¬ 
clusion;  Bq  =  const,  line  must  be  cut  first  for  our  pur¬ 
pose. 

(9)  B1  -  B3  , 

e3  =  hv*3e  ,  = 

Fig.  3-9  shows  e^  must  be  cut  first. .  This  is  interpreted, 
in  -  B^  terminology,  as  B^  =  const,  line  must  be  cut 

first  for  the  required  purpose. 

(10)  B2  -  B4 

The  same  analysis  requires  that  the  -  const,  line  should 
be  cut  first  as  the  equation  (3-5)  becomes 
F(  «n)  =  (  A4  *  B4)  e4  +  (  A2  -  B2  )  e2 

in  this  case. 

Seeing  these  rules,  it  is  desirable  to  classify  them  into  a  few 
groups  to  simplify  them,  as  given  in  the  following. 
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3.2  Summary  of  Enclosure  Criteria  -^line  cut  first  and  the  least  number  of 
cutting  -  points 


the  line 

the  line 

cut  first 

least  number  of 

cut  first 

least  number  of 

Case  I 

C  =  0 

cutting  points 

C  >  0 

.cutting  points 

B  -  B, 

B 

N  -  1 

B 

N 

o  1 

o 

0 

B  -  B„ 

B 

N  -  1 

B 

N 

o  3 

o 

o 

B1  ’  B2 

B2 

N  -  1 

B1 

N  +  1 

B2  '  B3 

B2 

N  -  1 

B3 

N  +  1 

B3  *  B4 

B4 

N  -  1 

B3 

N  +  1 

B4  -  B, 

B4 

N  -  1 

B1 

R  +  1 

These  six  pairs  permit  the  stability  check  in  its  strict  sense,  be¬ 
cause  they  have  no  singularities  at  £  =  0,  as  has  been  already  explained. 

The  next  group  consists  of  four  pairs,  all  even  -  even  and  odd-odd 
and  doesn't  have  the  stability  check  in  its  narrowest  sense,  because  all  of 
them  have  singularities  at  £  =  0. 


the  line 
cut  first 

least  number  of 

Case  II 

C  >  0 

cutting  points 

B  -  B- 

B 

N  +  1 

o  2 

Q 

B  B. 

B 

N  +  1 

o  -  .  4 

o 

B1  ’  B3 

B3 

N  +  1 

B2  '  B4 

B4 

N  +  1 

In  conclusion  it  can  be  said  that  the  line  for  6  -  const,  must 

even 

be  cut  first  when  stability  is  checked  in  any  case. 
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5*=  o,o 


Figure  3-1 
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A  =  12100 

o 

A,  =  6.4 

1  Figure  3-2 

a4=  1.) 
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W -  vo 


A  *  12100,  A,  =  682 
o  1 

A.  =  1.0 

4 


Figure  3-4 
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CHAPTER  IV 
SKETCHING  PROCEDURES 


4 . 1  Need  for  sketching  procedures. 

Mitrovic* s  method  is  a  very  convenient  method  as  faff  as  calcu¬ 
lations  are  concerned,  because  it  involves  the  use  of  real  numbers  only. 
Mitrovic  visualized  use  of  desh  calculators  with  secretarial  type  help, 
or  a  digital  computer,  or  an  analog  computer.  All  of  these  techniques 
take  time,  and  while  they  are  acceptable  for  detailed  studies,  the  en¬ 
gineer  making  preliminary  studies  prefers  a  technique  (such  as  root 
locus  or  Bode  diagram)  which  provides  reasonably  accurate  curves  with 
essentially  no  numerical  calculations  beyond  formulation  of  the  equations. 
Mitrovic* s  method  will  seldom  be  used  for  simple  problems  unless  such  a 
facility  is  available,  and  if  it  is  not  used  for  simple  problems  it  is 
not  likely  that  the  engineer  will  develop  the  familiarity  needed  to 
apply  it  to  more  complex  problems.  Fortunately  for  most  simple  cases  the 
Mitrovic  curves  for  £  -  0  can  be  sketched  rapidly,  and  for  £  4  0  only  a 
little  labor  is  involved.  Sketching  was  first  used  in  the  Elliott, 

Thaler,  Heseltine  paper,  but  was  not  developed  as  a  technique  for  gen¬ 
eral  use.  The  results  presented  in  this  chapter  are  intended  to  ex¬ 
tend  the  usefulness  of  sketching  methods. 

4.2  Basic  Manipulations  in  Sketching  The  Mitrovic  Curves. 

A  Mitrovic  curve  is  a  plot  of  two  equations  with  a  common  para¬ 
meter  which  is  the  frequency,  a)^  .  Since  the  two  equations  are  known,  and 
since  the  plot  (for  sketching  purposes)  covers  the  range  0  <  to  <  +  »  the 
locations  of  the  two  extreme  ends  of  the  curve  are  easily  found  by  letting 
=  0  in  both  equations  and  u>n  =  °>  in  both  equations .  With  the  two  ends 
located  evaluation  of  a  few  selected  points  on  the  curve  permits  sketching, 
and  fortunately  it  is  possible  to  locate  maximum  and  minimum  points  which 
indicate  the  extreme  excursions  of  the  curve.  Since  the  equations  for  the 
two  adjustable  coefficients  Bx  ,  B^  (also  the  coordinates  of  the  plot  1  ) 
are  known  it  is  easy  to  evaluate 


dB  dB 

=  0;  — 1  = 

da)  da) 

n  n 


0 


4  -  1 


These  ®qu&t ions  tan  be  evaluated  for  any  order  equation  and  for 
any  value  of  £,  but  for  sketching  purposes  £  =  0  Is  usually  wanted  since  it 
clearly  indicates  the  range  of  values  that  provide  a  stable  system.  Con¬ 
veniently  the  Mitrovic  equations  have  their  simplest  algebraic  fom  when 
£  =  0  so  that  numerical  formulation  of  the  equation  is  quite  elementary. 

Next  the  equations  must  be  -solved  for  their  roots  in  order  to  find  the 

values  of  a)  at  which  the  maxima  and  minima  occur.  For  characteristic 
tk 

equations  up  to  5th  order  (  and  in  some  cases  to  6th  order  )  the  poly¬ 
nomials  are  factorable  and  the  results  can  be  tabulated.  This  has  been 
done  for  a  reasonable  number  of  cases  and  the  results  are  tabulated  in 

Table  4-1.  Thus  values  of  to  for  the  maxima  and  minima  for  any  suitable* 

a 

combination  of  coefficients,  with  £  -  0S  can  be  found  by  referring  to 
Table  4-1  and  inserting  the  known  values  of  the  coefficient®  in  the  in¬ 
dicated  relationships.  It  is  then  necessary  to  return  to  the  original 
Mitrovic  equations  and  substitute  each  of  the  values  of  w  ,  thus  obtaining 
the  coordinates  of  the  maxima  and  minima  points.  The  Mitrovic  curve  for 
£  =  0  may  then  be  sketched  by  starting  at  ~  0  and  connecting  the  known 
points  in  order  of  increasing  w  ,  noting  that  each  known  point  is  also  known 

El 

to  be  either  a  maximum  or  a  minimum.  This  procedure  usually  provides  adequate 
information  for  preliminary  analysis  and  design  evaluation.  Where  a  slightly 
more  accurate  curve  is  desired  additional  points  may  be  calculated  in  the 
normal  fashion,  with  the  advantage  that  the  known  values  of  u)„  at  the  max¬ 
ima  and  minima  guide  the  choice  of  additional  to  values. 

u 

When  curves  are,  desired  for  £  4  0„  the  same  sketching  procedures 
apply i  and  Table  4-1  also  contains  a  number  of  entries  for  evaluation  of 
maxima  and  minims.  Note  that  when  £  4  0  the  Mitrovic  equations  are  more 
complex  algebraically  so  explicit  solutions  for  the  maxima,  and  minima  are 
restricted  to  equations  of  somewhat  lower  order. 

For  higher  order  equation® s  while  sketching  is  still  desirable, 
explicit  solutions  for  the  maxima  and  minima  are  not  available.  However, 
the  general  Mitrovic  equations  can  always  be  differentiated  and  the  diff¬ 
erentiated  forms  are  given  here  in  Table  4-2.  These  can  always  be  used  to 
find  the  maxima  and  minima.  Perhaps  a  simple  graphical  solution  is  fastest 
and  most  convenient?  the  proper  equations  are  selected  frem  Table  4-2  and 
into  each  selected  values  of  «  are  inserted,  obtaining  value®  for  a  B 

O,  &  *  Y 

Mr 
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Gordina tea 

o 

11 

3 

co  =  “ 

CO  or  CO  . 

.  max  rain 

BO  VS  Bl 

Highest 

Order 

all 

3 

4 

5 

6 

7 

3 

4 

5 

6 

Value 
of  £ 

0 

0 

0 

BO  , 

0 

+  00 

-  CO 

-  00 

+  ® 

+  ®! 

0 

0,  VA2/2A4 

O 

V 

fcl 

N3 

bo 

0,  L{A4  +  /a4-3A2A6  }/3Afe]^ 

Bl 

0 

+  °° 

•  oo 

—  00 

**  oo 

+  ® 

0 

0 

0,  Va3/2A5 

o,  Jkz/ik5 
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Continued 

co  or  co  , 

max  rain 

Highest 

Order 

7 

8 

Value 

of  t . . 

0 

BO 

— . 

0*  —  ^4  *  3A2A6  }^3A6 

% 

too 

complicated 

Bl 

°,  [{a5  +  Jk^  -  3A3A?  }/3A?  ^ 

!V 

too 

complicated 

TABLE  I  (continued) 


Coordinates 

10  =  0 

CO  =  ® 

CO  or  CO  . 

max  nan 

Bl  VS  B2t 

Highest 

Order 

all 

3 

4 

.  5 

6  . 

3 

4 

5 

6 

VaTue 

of  C  ^ 

0 

0 

0 

Bl 

0 

+  00 

+  00 

-  00 

0 

0 

0,  Va3/2A5 

0,  V  A3/2A5 

B2 

+  CO 

0 

V 

+  CO 

4-  GO 

i 

-  CO 

00 

<VVv4 

<aomV 

too  complicated 

Coordinates 

o 

H 

3 

CO  =  ® 

CO 

I 

or  CO  . 
nax  min 

B2  VS  B3 

Highest 

OEder 

all 

•  3 

4 

5 

6 

3 

4 

5 

6 

Value 

of  C 

0 

0 

0 

B2 

+ 

0 

4-  OD 

+ 

-  CO 

+  °° 

.  <W* 

<*A>k 

too  complicated 

B3 

+  ® 

0 

0 

+  ® 

+  00 

+  00 

+  “ 

(*1/A5)i 

(A1/A5)% 

TABLE  I  (continued) 


Coordinates 

o 

(1 

3 

03  = 

u)  or  to  . 

max  rain 

B3  VS  B4 

’  Highest 
Order  . 

all 

‘  4 

5 

r  * 

6" 

4 

5 

6 

Value 

of  C 

0 

0 

0 

B  3 

+  * 

0 

4.  00 

+  °° 

00 

(A1/  A5)^ 

(Aj^/  a5h 

B  4 

*  00 

0 

>-0 

+  ® 

^2  Ao/A2)  « 

•J  2  A  /A  ,  » 

►  O  £ 

too  complicated 

Coordinates 

O 

II 

3 

to  =  05 

to  or  w, 

max  min 

BO  VS  B3 

Highest 

Order 

all 

3 

4 

5 

6 

* 

T3 

4 

5 

6 

Value 
of  C 

0 

0 

0 

B  0 

0 

4*  00 

*  00 

•  00 

+  » 

0 

0  j  V A  2  /  2 

0,  7a2/2  A4 

too  complicated 

B  3 

4-  00 

0 

0 

+  « 

4.  00 

CO 

00 

(Al/A5)4t 

<a1/a5)^ 

coordinates 

CO  —  0 

CO  or  <0  . 

max  min 

Bl  VS  B4 

Highest 

Order 

all 

4 

5 

6 

4 

5 

6 

Value 

of  1 : 

,0 

0 

0 

B  1 

O 

4-  ® 

m  00 

-  CO 

0 

0,  VA3/2A5 

0,  y  A3/2A5 

B  4 

«■  ro 

0 

0 

■f  ® 

y  2a0/a2,  » 

^  2A0M2.  - 

too  complicated 

1 

* 


Coordinates 

o 

II 

3 

0)  “ 

w  or  w  . 

max  mm 

BO  VS  Bl 

Highest 

Order 

all 

3 

4 

3 

Value 
of  £ 

o<e<u 

0<C-S  0.5 

« 

o.5  <C^5  l 

0<C-5  0.5 

.  ’  .  . 

0.5  <C  <  0.71 

0.71^C^1 

h  * 

B  0 

0 

-  00 

-  00 

CO 

+  <» 

4*  « 

0,2A2/3A^P2(O 

B  1 

0 

+  ® 

*  00 

I  -  ® 

-  00 

4-  ® 

-a^2(0/a3<p3(C) 

Continued 

to  or  to  . 

max  min 

Highest 

Order 

4 

Values 
of  C 

B  0 

9,  {-3A3<P2  (C)  +/9A^<P2  (D  +  32A2A4<P3(0  }/8A4<P3  (O 

B  1 

{  -a3  <p3(C)  +  Va3  <p3(C)-3a2a4  <p2(D  <p4(C)  }  3A4  <p4CO 

N 


Coordinates 

e 

u 

o 

to  =  00 

to  or  '  to  . 

max  min 

Bl  VS  B2 

Highest 

Order 

all 

3 

4 

3 

4 

Values 
of  £ 

0  <C«i 

0  <C  -$1 

0  <£.<0.5 

0.5  <C<1 

B  1 

+  “ 

+  ® 

•*  CO 

»  00 

(Ao02(C)/2A3)1/3 

too  complicated 

B  2 

+  05 

+  ® 

+  CO 

—  CO 

(2Ao/A3^2<«1/3  . 

too  complicated 

Coordinates 

o 

11 

3 

CfJ  ~  = 

to  „  or  CO  - 
max  mm 

B2  VS  B3 

Highest 

Order, 

all 

3 

4 

3 

4 

Values 

of  C  < 

0<e  ^0.5 

0.5  <C«1 

■ 

0  <C 

0  <C^1 

B  2 

+  0£> 

*P  CO 

0 

+  00 

-2Aq  <P3(C)/Aj<P2(0,  “ 

too  complicated 

B  3 

m  CO 

*  CO 

0 

4*  oo 

3Aq  <P2  <C)/2A1(  » 

too  complicated 

% 


* 


Coordinates 

CO  =  0 

CO  -  03 

B3  VS  B4 

Highest 

Order 

all 

4 

5 

Values 
 of  £ 

0<£  ^0.5 

0.5  <£  <  0.71 

o.7i 

0< 

0<C<  1 

B  3 

•  CD 

-  CD 

+  ® 

0 

+  03 

B  4 

*  CO 

-f.  CD 

+  * 

0 

+  00 

Continued 

co  or  co  . 

max  min 

Highest 

Order 

4 

5 

Values 

of  g 

B  3 

{-2a1  <p3(C)  +  VaJp3<0-3aoa2  <p2(C)  <p3(C)}/a2  <p2(0>00 

too  complicated 

B  4 

{3AX  <P2(C)  ±  V9aJ  cP2(C)+32  AoA2  <P3(C)}/4A2<  « 

too  complicated 

JL 


Coordinates 

CO  =  0 

co  or  co  . 

max  min 

BO  VS  B3 

Highest 

Order 

all 

3 

'  4 

Values 

of  C  

0  <C  <0.5 

0.5  <C^1 

0  <£<0.5 

0.5  <£<1 

0  <£  <0.5 

0.5  <£  <0.71 

0.71  ^  1 

B  0 

0 

0 

+  00 

-  00 

•  00 

+  00 

+  00 

B  3 

+  00 

•  00 

-  00 

-  00 

+  00 

•  00 

+  00 

Continued 

oo  or  03  . 

max  min 

Highest 

Order 

3 

4 

Values 
of  £ 

B  0 

Ax  (p2(D/ 2  A2  . 

too  complicated 

B  3, 

2  A1/A2  <P2(0,“ 

too  complicated 

Coordinates 

CO  =  0 

U}=  « 

CO  CO  . 

max  or  min 

Bl  VS  B4 

Highest 

Order 

all 

4 

4 

Values 

Of  £ 

0<£  <0.5 

0.5  <C  <  0.71 

0.71  <C  <1 

0  <e  <  0.5 

0.5  <C<1 

B  l 

4*  ® 

<m  OP 

•  00 

too  complicated 

B  4 

—  00 

+  » 

4,  00 

0 

0 

too  complicated 

TABLE  I  (continued) 


Coordinates 

Ui  —  0 

GO  =  ® 

BO  VS  B2 

Highest 
Order  

all 

3 

4 

Values 
of  C 

o  <C^i 

0  <C  <0.5 

o.s^c^i 

0  <C<0.71 

0.71  <£.<1 

B  0 

0 

•  00 

-  00 

4-  oo 

+  00 

B  2 

00 

-  00 

+  00 

+  ® 

-  00 

> 


* 


Continued 

co  or 

max 

U)  . 

min 

Highest 

Order 

3 

4 

Value 
of  C 

BO 

h1cp2  <0/2A2 

too  complicated 

B2 

2  Ax/A2  <p2  (O,  00 

too  complicated 

TABLE  I  (continued) 


Coordinates 

oj  =  0 

03  =  ® 

Bl  VS  B3 

Highest 

Order 

all 

3 

4 

5 

Values 
of  C 

0  <C  <0.5 

0 . 5 .<C 

0  <C^i 

0  <C  <0-5 

0.5^C$1 

0  <g£0.71 

0.71  <C  <1 

Bl 

-  oo 

■f  03 

+  * 

-  00 

-  CD 

4.  00 

+  1:0 

B3 

m  00 

—  CD 

0 

*  00 

+ 

4  00 

•  00 

* 


Continued 

LC  or  CO  . 

max  min 

to  or  oj  . 

max  min 

Highest 

Order 

3 

4 

5 

Values 
Of  c  . 

B  1 

too  complicated 

7a0(p3«>/*2 

L{a2±7a^+  12  A0<P3(C)}/6aJ1/Z 

B  3 

.  f/3A  /A_,  <” 

0  i 

[{-A2  t  ■hi  -  12AoA3  V3(C)}/2A4]l/2 

too  complicated  j 

lAJldi  4*  \  4  NU^U  J 


Coordinates 

u)  =  0 

CO  =  » 

oo  or  to  . 

max  min 

B2  VS  B4 

Highest 

Order 

all 

4 

4 

Values 

Of  c 

0<C<0.71  0.71  <£< 

* 

0<C«1 

B2 

+  OCJ  -  00 

+  00 

too  complicated 

B4 

+  00  +00 

0 

V 


* 


Coordinates 

to  -  0 

cu  =  » 

B0  VS  B4 

Highest 

Order 

all 

4 

J 

5 

Values  ' 

of  c 

0  <£<0.71 

0.71  <0^1 

0  <£<0.71 

0.71  <£<l 

0  <£<0.31 

O.31^C<0,71 

0.71  <C<1 

B0 

0 

! 

0 

+  ® 

-  00 

*  00 

*  OO 

■+  00 

B4 

-  00 

+*> 

0 

0 

-  00 

+  “ 

+  * 

taju£  i  ^continued; 


Continued 

%in 

Highest 

Order 

4 

5 

Values 

of  c 

BO 

■ao<p2(C)  ±Va2<p2  <C)  +  SA^CC)  }/3A3 

too  complicated 

B4 

\  -a2<p2(D  ±7a2  (C)  +  3  Aj_  a^3(o}a3<p3(C) 

too  complicated 

at  selected  cu  values,  A  plot  is  then  made  of 
n 

Each  curve  crosses  the  to  axis  at  values  of  to 

n  n 

and  minima. 


■  vs  lo^  and  vs 

corresponding  to  the  maxima 


4*3  Tabulation  of  Results  and  Suggested  Procedures. 

(Tables  containing  useful  relationships  are  collected  in  the  pages 
following. 


Suggested  Procedures:  (these  are  largely  a  summary  of  the  discussion 
in  section  4*2) 


1.  From  the  characteristic  equation  determine  the  desired  Mitrovic 
equations  for  £  =  0. 

2.  Substitute  in  these  equations  GO  -  0  and  to  -  ^  to  locate  the 

-  n  n 

ends  of  the  curve* 

3*  Inspect  the  Mitrovic  equations  to  see  if  a  ready  solution  exists 


for  zero  points;  ie*,  for  what  values  of  to  are  B  =  0,  B  =  0? 

n  x  y 

dBx  =  0, 

<3uT  dto 


4.  Obtain  the  equations 


=  0  and  solve  for  values 


of  co 


n 


n 


5*  Substitute  the  values  obtained  in  3  and  4  into  the  Mitrovic 

equations  to  locate  the  B  ,  B  coordinates. 

x  y 

6.  Plot  the  known  points  and  sketch  the  curve. 

4.4  Illustrations  of  the  Sketching  Techniques 
Example  I 

For  a  third  order  equation  vs  £  =  0 
F(s)  =  s3  +  2  s2  +  B  s  +  Bo 


C  =  0 


B  =  2  uf 
o 

B1  =  CO2 


(4-1) 


From  the  relationships  in  the  table,  Bq  and  are  both  equal  to  zero  at 

O)  =  0;  they  have  no  maximum  or  minimum  points;  both  end  in  +  ®.  Thus  for 

equation  4-1,  a  sketch  of  B  -  B,  curve  is  shown  in  Fig.  4-1.  Notice: 

o  i  > 

all  of  the  roots  of  F(s)  have  a  damping  ratio  £=0.5  if  B  =  1,  B  =  2. 

0  1 
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TABLE  4-2 


DIFFERENTIATED  FORMS  OF  MITROVIC9  EQUATIONS 
FOR  USE  IN  EVALUATING  MAXIMA  AND  MINIMA 


BO  VS  B2 

=  {-2  A2u+  3  A,*2  «)  J  +  ...  +  n  An  Vl(D  w"-1] 

=  *2*2  <£)  +  2  *3^3 (C)  «+•••+  <“-D  *„«>„  «>  “P'2 

BI  VS  B2 

-  -A„  »>2  (O  5-  +  2  Aj  u  -3  A4.«>2  (C)  m2  -  ...  -(n-l) 

An V2  <0  “'2 

W  ”  *2  A3,)2  (D  +  2  A4*3  <£>  “+  •••+<"-2> 

(C)  A)  n‘3 

B2  VS  B3 

d  co  =  ”2  Ao  ^3  ^  ^  ”A1  ^2  +  2  A4  w  +  •  •  •  +  (n-2) 

An  *n-3  «>  W"'3 

=  3  Ao  ^2  “3*  “2  Al  “3  +  A4  <P2  (O  +  . . .  +  (n-3) 

n-4  W 
A«  ^n-9  W 

n  n-z 

B3  VS  B4 

=  _3  Ao  *4<«-j*  _2Ai  <p3  (o~j  -A2  ^2  (C)4r2  +  2A5U> 

-<n  ’  3>  A„<Pn-4(°  ^ 

~f^T  =  4  Ao  <P3  (C)"^5-  +  3  A1  <P2  (0”^7T  “2  A 2~^  +  A5  <P2  (C)  + 

...  +  (n-4)  Ancpn_3  (O  C0n'5 
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BO  VS  B3 
d  B  ( 


TABLE  4-2  (Continued) 


d  00 

*3<C>  L 

d  B  3 

-  1  f 

d  oo 

'  V3(C)  1 

VS  B4 

d  B  1 

i 

d  oo 

'  <p3(0 

d  B  4 

1 

d  oo 

"  <P3( C) 

VS  B2 

d  B  0 

d  OO 

"  <P2(C)  [ 

d  B  2 

-1 

d  oo 

'  <P2(C) 

VS  B3 

d  B  1 

=  1 

d  oo 

<P2  (D 

d  B  3 

l 

d  60 

"  ^(C) 

(D  con_1  ] 


<C)  “"'4  ] 


.  +  (n-1) 


W4  «>  "‘2  ] 

T  4Ao  +  2  “T  +  A,  <P2  2  ”  •  •  • 

L  U  w  w 

An  Vl  ^  1 


-  (n-4) 


(?)  oo11'1  ] 


CO2 


(n-2) 


An^n  «>  “~3] 


d  B  1  =  1  [  A  <p3  (C)  +  A2  -3  A4  w2  +  4  A  <p2  (O  <J  +  . . .  + 

d  CO  <p„  (O  00 

<~l>  *„  V3«>  “”'2 1 

4  A  A„ 


4  “  2  "A4  ^3  (C)  “  2  A^  <p^  (()  oo  -  • . .  "(n-3) 

c^.r<o  ‘-n-4  ] 
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B2  VS  B4 


TABLE  4-2  (Continued) 


d  B  2  

1 

d  (0 

<P2(C) 

d  B  4  

1 

d  CO 

ro 

VS  B4 

d  B  0  

1 

d  co 

<P4(C) 

d  B  4  

1 

d  W  <P4(C) 


[  2  Aq  <p4  (O  -j-  +  Ax  <p3  (O  +  . . .  +  (n-2) 

An  *>n-4  «>  "'3  ] 

[-4Ao^2«)-^+  ~3T  -J---(n-4) 

A„  ”„-2  «>  “"'5  ] 

[  Ax  <P3  (O  +  2  A2  <P2  (C)  to  -3  A3  CO2  +  5  A5  co4 

*  •••  •  "*«*.•*  «>  u”'1  ] 

[  -  1T+  2  A2*2  «>  ^  +  A3*3  «>  ~J  -*S*,S«> 

...  -  (n-4)  An  «3n  (O  w*‘5  ] 
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Example  II 


F(s)  = 
s  =  -1, 


3  2 

s'1  +  s  +2s+l= 

.1  +  j  JT 
2 


(s  +  l)(s2  +  s  +  I) 


For  a  fourth  order  equation  B  vs  B. i  ^  =  0  and  £  =  0.2. 

4  3  °  2  L 

F(s)  =s  +4. 2s  +6. 6s  +B^s+Bq 


For  C  “  0, 


B  =  A,  u>2  -  A,  w 
o  2  4 


B1  =  A3  40 


(4-2) 


The  table  indicates  thatboth  B  and  B.  begin  from  zero;  B  ends  at 

O  X  ~p — - — o 

whereas  B-^  ends  at  +  Bq  has  a  maximum  at  0)  =  7a2/2A  .  Also  equation 

(4-2)  indicates  that  B  =  0  at  to  =  JhTfh.  .  Thus  : 
v  o  o  Z  4 


to  for  B  -  76.6/2  =  1.8 
max  o 


to  =76.6  =  2.6 
o 

Bo(1.8)  =  10.9  B1(1.8)  =  13.9 

Bq(2.6)  =0  Bx  (2.6)  =  27.7 

The  curves  are  shown  in  Fig.  4-2. 

For  £  =  0.2, 

bq  =  a2  to2  -  a3  <p2(C)  to3  -  a4  <p3  (C)  to4 

2  3  (4-3) 

B1  =  A2  ^2^  u  +  A3  ^3^  W  + 

The  table’s  indications  are:  both  Bo  and  B^  are  equal  to  zero  at  to  =  0; 
both  Bq  and  B^  are  equal  to  -  ®  at  to  =  <®,  if  £  =  0.2;  each  curve  may  have 
a  maximum.  The  numerical  values  of  the  <p (£)  functions  are; 
ip2(0.2)  =  0.4 

(p3(0.2)  =  0.84 

tp4(0.2)  =  -0.736 

The  positive  values  of  comav  in  the  table  are  the  only  ones  which  can  make 
Bq  and  B^  maximum,  and  they  are: 
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(0  for  B 
max  o 


=  -3  x  4.2  x  0.4  +  7(3  x 

=  1.4 


(3  x  4.2  x  0.4)  +  32  x  6.6  x  o.84  = 

8  x  0.84 


to  for  B. 
max  1 


=  -4.2  x  0.84  -  V( 4.2  x  Q.84)2  +  3  x  6.6  x  0.4  x  0.736  = 


(4.2  x  0.84)' 

-3  x  0.74 


=  3.5 


The  corresponding  curves  are  shown  in  Fig.  4-3.  Notice:  all  of  the  roots 
of  F(s)  have  a  damping  ratio  £j  =  0.6,  if  Bq  =  2,  Bj  =  5.4.  This  point 
(5.4,2)  is  shown,  in  Bq  -  curves  for  £  =  0  and  £  =  0.2,  to  be  in  the 

enclosure. 

Example  III 

For  a  fifth  order  equation,  Bq  vs  B^,  £  =  0 

F(s)  =  s5  +  7  +  18  s2  +  23  s2  +  B.  s  +  B 

i  o 

2  4 

B  =  AjU  -  A, CO 
O  Z  H 

2  4 

Bi  ”  A^kJ  m  A^oj 


(4-4) 


According  to  the  indications  of  the  table,  both  Bq  and  B^  are  zero  at  to  =  Os 


B  =  B.  = 
o  1 


at  to  =  each  has  a  maximum  at 


Wmax  (for  V  =  VA2/2A4  =  ^23/14  =1.3, 


‘Hnax  (f°r  Bl}  *  Va3/2A5  =  Vl8/2  =  3. 

Equation  4-4  shows  that  B  is  zero  at  to  =  -/A„/A,  =  1.8;  B.  is  zero  at  to  = 
7-3 -  o  o  2  4  1  O 

=  vAJ/A_  =  4.2,  Thus: 
o 


Bo  (1.3)  =  18.9 
Bq  (3)  =  -360 


Bl  (1.3)  =  26.8 


Bo  (1.8)  =  0 

,  B  (4,2)  =  -1854 
o 


Bj  (3)  =  81 
Bx  (1.8)  =  48,3 
Bl  (4.2  =  0 

The  curves  are  shown  in  Fig.  4-4.  Notice:  All  of  the  roots  of 
F(s)  have  a  damping  ratio  £  =  0.5  if 
Bq=  6 

B1  =  17 
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This  point  (17,6)  is  shown  in  Fig.  4-4  to  be  in  the  stable  domain. 


Example  IV 

For  a  seventh  order  equations  8q  vs  B^,  £  =  0 

F(s)  =  s  +  9  s  +  33  s3  +  66  s  +  81  s  +  63  s  +  BjS  +  Bo 

2  4  6 

B  =  A  cu  -  A,uj  +  A, co 
o  2  4  6 

2  4  6  <4'5> 

B1  =  Al"  *  Aj“  +  a7“ 


Table  4.1  shows  that  both  B  and  B.  are  zero  at  to  =  Oj  B  =  B,  =  ®  at  oj=  ®; 

o  1  o  i 

both  can  have  one  maximum  and  one  minimum 9  but  are  not  required  to  have  such 
singular  points  other  than  at  co  =  0. 

Table  4.2  shows: 


d  B  ,  - 

— =  2A„o)  -  4A,  oj  +  6A,l<j 
dw  2  4  6 

(4-6) 

^  35 

“53"  =  ^3^  "  4A5w  + 


For  this  case  B^9  has  three  real  roots; 

U>  =  0,  <2  =  132  ±Vq32)2  -  4  x  27  x  63  =  4.2  and  0.35 

54 

W  =  0,  0.74,  and  2.05.  .B*  also  has  three  real  roots;  uj =  0,  1.15,  and  4.5. 


The  curves  may  have  a  shape  as  in  Fig.  4-5.  From  these,  the  shape 
of  Bq  *  is  apparent  from  the  sequence  of  frequency  for  the  maximum,  min¬ 
imum  and  zero  points  as  shown  in  Fig.  4-6.  Note  that  either  (a)  or  (b)  sketches 
are  possible,  because  the  magnitudes  of  £q  and  have  not  been  calculated. 

If  the  critical  values  of  io  are  substituted  into  the  Bq  and  B^  equations,  the 
coordinates  of  the  points  iare  determined  and  the  sketch  can  be  made  accu¬ 
rate.  Thus 


£0 

0.74 

2.05 

1.15 

4.5 

2 

GO 

0.55 

4.20 

1.33 

20.7 

4 

to 

0.30 

17.65 

1.77 

428 

6 

w 

0.16 

74.1 

2.35 

8870 

4  -  21 


and  Che  coordinates  of  the  critical  points  are; 


Bq(0.74)  =  172 

Bo(1.15)  =  -15 

Bo(2.05)  =  -237 

B  (4.5)  =  52905 
o 


Bj (0. 74)  =  35 
B1(1.15)  =  48 
B1(2.05)  =  -170 
B^(4.5)  =  -3553 


Thus,  the  B  vs  B.  curve  is  shown  to  be  as  in  Fig.  4-7.  Notice:  All  of 
°  > 

the  roots  of  F(s)  have  a  damping  ratio  £  =  0.5,  if  B^,  =  29,  Bq  =  6. 

F(s)  =  s7  +  9  s6  +  33  s5  +  66  s4  +  81  s3  +  63  s2  +  29  s  +  6 

=  (s  +  l)2(s  +  )2  (s  +  3)(s2  +  s  +  1). 

This  point  (29.6)  is  shown  to  be  inside  the  stable  domain. 

Example  V 

For  the  fourth  order  equation,  B2  vs  Bj,  £  =  0. 

F(s)  =  s4  +  B3  s2  +  B2  s2  +  5.4  +  2 

B  =  A  -Kr  +A  W 
2  o  2  4 

CO  H 

b3  •  h-h 

CO 


(4-6) 


The  indications  of  the  table  are  that  both  B2  and  B^  are  »  at  co  =  0J  B2  =  ® 
B^  =  0  at  a)=  ®;  B2  has  a  maximum  at 

C0=  (Aq/A4)1/4  =  (2)l/4  =  1.18 

Thus: 

B2  (1.18)  =  2.82, 

B3  (1.18)  =  3.83 

and  the  curves  are  shown  in  Figs.  4-8  and  4-9.  Notice  all  of  the  roots  of 
F(s)  have  a  damping  ratio  £  =  0.6  is  B2  =  6.6,  B^  =  4.2 

F(s)  =  s4  +  4.2  s3  +  6.6  s2  +  5.48  +  2 

=  (s  +  l)(s  +  2)(s2  +  1.2  s  +  1). 

The  point  (4.2,  6.6)  is  shown  in  the  graph  to  be  lying  in  the  stable  domain 
Example  VI 

For  a  fifth  order  equation,  B^  vs  B2,  £  =  0 
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The  table  shows  that  =  0,  B2  =  +  ®  at  to  =  0;  Bj  =  -®,  B2  =  +  00  at  to  = 
B^  has  a  maximum,  and  B2  has  a  minimum  at: 

"max  =  •/I7*rS  =  /I57r=  3 

"min  *  <W1M  '  <6/7>1/4  =  °'96 

Equation  4-7  shows  that  B^  is  zero  at 

to  =  s/Ajk'  =  4.24 

o  3  5 

Thus,  the  critical  points  are  at: 


B1(0.96)  =  15.7  B2  (0.96)  =  12.9 

B^ (3)  =  63.6  B2  (3)  =  63.6 

B1(4.24)  =0  B2  (4.24)  =  126.3 

The  curves  are  shown  in  Fig.  4-10 
Example  VII 

For  a  fourth  order  equation,  Bj  vs  B2,  £  =  0.6 
4  3  2 

F(s)  =  s  +  4.6  s  B2s  +  BjS  +  2 

B1  *  Ao<P2(°4r  +  As"2  •  A4  *^2  <C)  "3 

B2  =  A0  ij  +  Aj  <P2(t)  CO  +  A^  <P3(C)  C02 


Table  4-2  gives  the  following  equations,  if  the  corresponding  values 
of  and  (C)  are  substituted; 

d  B, 


-57-^2  <2'4  '  “3  +  3-“»4) 

CO 

^B2  1  34 

=  -  -j-  (4-5.52  uT  +  0.88  to*) 

or 


(4-9) 


From  the  Table  4-1,  it  is  seen  that  B^  =  B2  =  +  ®  at  go  =  0; 


=  B2  =  -®  at  to  -  ®.  When  £  =  0*6,  the  maximum  or  minimum  frequency  is 
not  available  through  the  routine  methods*  Even  their  existence  is  not 
readily  checked* 

Equation  4-9  can  be  used  to  improve  the  situation*  It  is  seen 
from  figure  4-11  that  there  may  be  two  points  where  the  difference  of  the 
curves  is  just  2*4,  since  by  the  differential  calculus  it  is  easily  shown 
that  the  difference  between  the  two  curves  attains  its  maximum  value  of 
16.3  at  to  =  1.91.  So  there  must  be  one  frequency  before  1.91,  and  an¬ 
other  after  1.91,  which  realize  these  conditions.  It  is  rather  easy  to 
find  these  values  by  trial  and  error,  and  they  are  found  to  be $ 
to  “  0.7  and  2.5 

By  reasoning  of  the  whole  figure,  the  smaller  one  must  be  the  minimum 
point,  and  the  larger  one  the  maximum. 

Through  the  same  analysis  applied  to  the  other  equation  of 
4-9,  it  is  seen  that  to  =  0.95  and  6.2  are  the  required  ones.  Thus,  this 
second  relationship  provides  Fig.  4-12.  and  B2  are  then  plotted  as  in 

Fig.  4-13  and  the  Mitrovic  curve  is  sketched  in  Fig.  4-14 

If  the  critical  values  of  to  are  substituted  into  equation  (8) , 
the  graph  becomes  more  accurate | 


Bl(0.7)  =  5.21 

B1(0.95)  =  5.51 
Bx(2.5)  =  11.22 
B1(6.2)  =  -109 

These  are  deduced  by  using  the  table g 


B2(0.7>  =  7.65 
B2(0.95)  =  7.05 
B2(2.5)  =  11.3 
B2(6.2)  =  17.4 


00 

0.7 

0.95" 

2.5 

6.2 

J  1 

0.49 

0.90 

6.3 

38.4 

to3 

0.34 

0.86 

15.6 

238 

1/co 

1.4 

1.0 

0.4 

0.161 

1/to2 

2 

1.1 

0.157 

0.026 
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A  more  accurate  sketch  of  the  vs  B2  curve  is  then  available,  as  shown 
in  Fig.  4-15.  Notice:  All  of  the  roots  of  F(s)  have  a  damping  ratio 
C  =  0.8,  if  Bl  =  6.2,  B2  =  7.8. 

4  3  2 

F(s)  =  s  +  4.6  s  +  7.8  s  +  6.2  s  +  2 


=  (s  +  l)(s  +  )(s2  +  1 
This  pair  is  shown  to  be  in  the  enclosure. 


6  s  +  1) . 
as  should  be. 
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Figure  4-13 


Figure  4-14 
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CHAPTER  V 

THE  DIGITAL  COMPUTER  PROGRAM 

The  aims  of  the  present  programming  are; 

(i)  To  print  out  the  calculated  values  of  any  pair  of  Bg  and 

B  : 
m 

(ii)  To  make  a  graphical  plot  of  this  pair  on  section  papers  for 
each  value  of  £. 

The  first  is  done  by  the  usual  digital  machine,  and  the  second  is 
done  by  the  digital  analogue  converter  attached  to  the  system. 

Some  safe-guard  statements  had  to  be  included  in  the  programming 
owing  to  the  generalization  of  applicability  to  airy  pair,  in  addition  to  the 
original  purposes  of  getting  graphs  and  calculated  results  for  varying  £. 

These  will  be  shown  in  the  next  section. 

The  practice  of  programming; 

What  is  needed  to  run  the  machine  is  just  to  punch  the  cards  according 
to  the  indication  shown  in  the  beginning  part  of  the  programming; 

1)  User's  name  and  program  identification; 

2)  Output  specification. 

Blank  card  is  for  the  case  both  print  and  graph  are  re¬ 
quired  . 

Punch  "graph1.1  if  graph  only  is  required. 

Punch  "print"  if  print  only  is  required. 

3)  Order  of  equation. 

4)  Coefficients  of  equation,  Aq  . . . ,  Ar,  including  the  Ae  and 

A  which  are  assumed  to  be  variable, 
m 

5)  Values  of  e  and  m,  such  as  1  and  2  for  B^  -  B£  case. 

The  values  of  Bg  and  B^  vary  in  two  ways;  one,  through  the 

changing  value  of  £,  and  another,  through  a). 

Usually  there  are  no  difficulties  in  changing  the  values  of  £ 
from  0  to  1  by  an  arbitrary  chosen  step  if  the  difference  of  e  and  m  is 
not  even.  If  it  is  even,  B  values  become  infinite  and  the  machine  stops, 
as  £  takes  the  value  zero.  The  same  situation  occurs  if  the  difference  is 
a  multiple  of  three  when  £  is  equal  to  0.5.  These  difficulties  are  auto- 
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raatically  avoided  by  making  the  machine  choose  the  adequate  starting  values 
and  step  width  of  £  in  the  programming ■ 

The  problem  of  how  to  find  out  the  proper  upper  limit  of  10  comes 
next.  It  seems  very  logical  to  set  the  limit  around  the  value  of  u)  which 
makes  the  last  term  of  each  B  dominate,  or  become  larger  than  the  sum  of 
other  terms,  since  this  actually  assures  that  the  curve  is  plotted 

into  the  last  quadrant  where  it  ends  at  co  -  infinity.  This  is  also  necessary 
from  the  practical  points  It  sometimes  occurs  io  certain  pairs  of  and 
Bm  with  certain  values  of  £  that  the  curve  re- appears  in  the  first  quadrant 
at  the  larger  frequencies  and  divides  the  enclosure  or  stable  domain,  and 
changes  a  part  of  it  to  lie  out  of  enclosure  or  unstable.  Rather  a  lengthy 
part  is  spent  in  the  programming  to  secure  the  effective  results.  But  in 
many  cases,  the  plotted  graph  of  B&  -  becomes  inaccurate  by  extending 
the  value  of  to  this  much,  though  the  graph  scaling  is  automatically  set  for 
each  value  of  £  in  order  to  get  a  full  resolution,  since  the  values  of  them 
at  the  start  and  at  the  end  differ  tremendously.  On  this  occasion,  it  is 
rather  recommended  to  try  to  set  the  upper  limit  of  to  to  just  half  of  the 
value  indicated  by  the  machine.  This  decreases  the  values  of  the  highest 
B^  and  to  nearly  one-sixteenth  in  the  case  of  fourth  order  equations. 

It  is  a  good  practice  to  decrease  the  values  of  to  after  getting  the  printed 
results  to  get  graphs  with  better  resolution  whenever  it  is  necessary.  It 
also  improves  the  starting  part  of  the  printed  results,  and  thus  sometimes 
shows  the  maximum  or  minimum  at  the  very  lower  frequencies  hidden  in  the 
original  setting.  Anyhow,  it  is  necessary  to  know  the  tendency  of  the 
curve  at  o)=  0  and  at  a)  =  “  for  correct  understanding  of  the  results,  and 
it  is  very  easy  to  get  these  knowledges  by  a  glance  at  the  -  B  equa¬ 
tions  a 

To  improve  the  printed  results  is  easily  done  by  omitting  certain 
statements  which  are  inserted  to  make  the  machine  print  out  only  every  tenth 
out  of  the  900  calculated  results  in  th£  memory  bo:*  required  for  the  graph. 

It  may  be  said  that  many  practical  applications  come  out  of  this 
programming,  and  if  some  applications  are  beyond  the  present  programming, 
it  may  be  easily  attained  by  just  manipulating  a  few  statements s  For  an 
example,  if  the  behavior  of  the  curve  of  -  Bq  at  (  -  0  is  needed,  it  is 
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easily  attained  by  making  the  starting  value  of  £  negative  and  the  ending  value 
positive  just  around  the  origin  and  thus  getting  the  full  account  of  the  curve 
tendency  there.  Actually  just  two  statements,  namely,  statement  74  and  49  are 
to  be  manipulated  in  this  case. 
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CHAPTER  VI 


A  PROCEDURE  FOR  COMPENSATION  DESIGN 


6-1  Cascade  Compensation 

Throughout  the  procedure  of  compensation  design  given  heres 

K&,  or  Kp  is  assumed  to  be  the  same  as  the  original  system. 

It  is  assumed  that  a  system  with  a  transfer  function 
m 


G  (s)  = 


Ko  *«•  +  *1> 


»N  jH.'.s  +  pt) 


TT 


(6-1) 


is  to  be  compensated  by  a  filter  with 

0  (S)  =  Kg  5  +  * 
c  S  +  p 

to  have  £  >  ^  *  Then  the  characteristic  equation  is 


s,,+n+i+*N+n<p>sN+n+ 


o  0  4*  A 


(6-2) 


,m+1(P,  Kc)sm+  V 


+  Aj_  (P,  Kc)  s  +  Aq  (P,  Kc)  =  0 


(6-3) 


Actually  the  last  term  is 


A  (P,  K  )  -  K  .  P 
o  c  o 

as  K  *s  are  constant* 
v 

It  is  very  remarkable  that  any  of  A^  (P)  and  A^  (Ps  Kc)  is  a 
linear  function  of  P  and  P  and  (This  can  be  seen  easily*)  This  im¬ 

plies  that  the  gradient  of  B  or  B  has  some  value  independent  of  the  value 

e  m 

of  P,  determined  only  by  values  of  £  and  and  Bc  and  are  both  linearly 
expressed  in  A^  and  A^.  This  is  seen  as  follows: 


Be  <P'V  -1,5,11^1  <p-  V 


Bm  <P'Kc>  «>  *rA=  <P>  Kc> 


(6-4) 


The  effect  of  the  change  of  P  and  K  values,  fixing  £  and  00  is  shown  as* 

c  Ifc 


1  f  d  it  ®  K  ~\ 

A  B  =  £  <p,  (O  ur1  -r— ^  AP  +  ~^AK 

e  i.i.k  1  nLop  d  K  cJ 

3  J  1  C 


9.  A 


9.  A 


(6-5) 


AB  =  £'  (p  (C)  ojb  r  ~  a  P  +  zrir  ak  1  , 

m  a,b,c  a  ia  L.  9_  d  c  J 
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So  and  A£  are  all  linear  in  P  aod  K 


9  it  is  seen  that  both  A  B  ,  and 

C  6 


A  B  depend  only  on  the  differences  of  p  and  K  a 


m 


it  is  seen 

that 

Sap 

9  A 

A  k  = 

e  op 

A?  +  hT 

C 

AKC 

(6-6) 

a  a 

9a 

*  *  m 

A  A  -  5— 

APK-S 

A  K 

m  S  p 

9-k 

C 

c 

If  an  arbitrarily  chosen  P  and  make  the  M  (A^*  A^)  point  lie 
on  or  inside  the  Hitrovic  curve  B  -  B  with  specified  £  -  £  »  a  desired 

£  ID  O 

damping  has  been  obtained  and  the  design  procedure  ends  *  But  this  seldom 

occurs  in  involved  problems.  In  such  cases  P  and  must  be  varied  so  as 

to  make  a  new  M?  (A\  As)  point  lie  on  or  inside  the  new  B*  *  Bfl  curve . 

e  m  e  m 

This  is  attained  by  makings  at  a  certain  value  of  a)  ,  sav  at  lc  -  oj  f  M" 

*  n  o  o 

locus  cut  “  B11  curve.  This  is  mathematically  shown  as 
e  m 


A'  =  B* 
e  e 

A*  =  B* 


m 


m 
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where 


Al  -  A _  + 


A1  =  A  + 
m  m 


B9  =  B„  + 


B9  =  B  + 
m  m 


9 

A 

e 

9 

P 

9 

A 

m 

a- 

P 

9- 

Be 

9 

P 

9 

Bm 

9 

P 

API 

A  P  + 
A  P  + 


*  Ae 

Tk 


9  A 


m 


m  A  F  + 


9 

K 

9- 

Be 

9 

K 

c 

9 

B 

m 

AK 


A  K 


A  K 


(6-8) 


9  K 


All  A  ,  A  ,  B  and  B  are  previously  calculated  numbers.  Trial  values 
e  m  e  m 


of  P  and  determine  them.  Substituting  (8)  into  (7) * 


m 


,  9  B 

9  A  N 

.  9 

B 

a 

A 

A  P  ( 

.ITT 

e  ^ 
a  p  y 

1  +  A  Kc  ( 

„  5" 

e 

Kc 

*  9 

e 

kC 

^  9  B 

9  A  N 

,  9 

B 

‘  9 

A 

A  P  ( 

m 

m  ^ 

t  +  A  K„  ( 

m 

e 

v  9  p 

”  9  p  y 

s.  9 

K 

c 

"  9" 

K 

c 

) 

) 


(6-9) 


From  this  A  P  and  A  Kc  are  given  by 
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(6-10) 


AP  = 

D 


A  K  = 
c 


Kc 


where 


-(  a 


( 


9  p 
a  b 


a  a 

_ e 

3  p 

9  K_ 


-n  >■  a  Bm 

X  " 


a  a 


9  K, 


r)  -( 


a  b  a  a 


_ m 

a  p 


a  p 


S') 

p  J 


) 


9  b  9  A 


m 


11 

a* 

Q 

e  V 

-  9  K 

c 

dkc  = 

Cm 

/»*« 
V  9  p 

c  = 

e 

A 

e 

-  B 

e 

C  = 
m 

A 

m 

"  Bm  ' 
m 

r)  -<*( 


9  B  9  A 


a  k 


a  k 


) 


a  a  ^  -9b 

e)-<X  " 


9  P 


9  A 

m 

a  p  a  p 


) 


It  must  be  remembered  that  all  of  the  terms  in,  arty  bracket  are  numbers 
as  £  and  LO^  are  specified. 

Finally,  if  the  newly  found  P'  and 

P'  =  p  +  A  P 

K  ’  =  K  +  A  K 

c  c  c 

are  physically  acceptable,  namely  both  P”  and  are  positive  and  the  mag¬ 
nitude  of  Kc  is  reasonable,  then  the  design  is  completed  and  z  is  found  to 
be  z  =  P'/k£.  If  the  solution  is  not  acceptable  this  simply  implies  that  a 
multiple  section  filter  must  be  used. 

If  this  happens,  one  tries  to  find  values  of  P  and  Kr  as  favorable 

c 

as  possible  for  the  first  stage*  The  next  stage  is  just  a  repetition  of 
the  procedure  given  here*  This  process  can  fee  carried  on  until  the  spec¬ 
ifications  are  met* 

It  also  must  be  noticed  that  although  the  above  mentioned  method 
is  very  general,  it  is  sometimes  found  that  a  simpler  treatment  is  applicable 
and  preferable. 

If  the  order  of  the  whole  compensated  system  is  high  so  as  to  have 
more  than  one  pair  of  complex  conjugate  roots  whose  real  parts  are  in  the 
same  order  of  magnitude,  the  above  mentioned  procedure  must  be  followed 
rather  strictly* 
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But  if  it  is  known  that  there  i®  just  one  pair  of  complex  con¬ 
jugate  roots,  or  of  dominant  complex  conjugate  roots,  then  more  simplified 
procedure,  is  applicable. 

In  this  occasion  it  is  justified  to  locate  M  Affi)  point  on 

any  part  of  the  Bg  -  Bffi  curve  in  the  first  quadrant.  It  is  not  necessary 
to  place  them  on  or  inside  the  enclosure  part  of  the  curve.  Then  the 
problem  degenerates  to  just  finding  out  the  values  of  P  and  which  makes 


A  =  B 
e  e 

A  =  B 
n  m 


\  «p-  K^)  takes  the  following  forms 

\  CP.  K0>  ■»kPtfkle  +  V 

where  and  are  all  given  real  numbers. 


(6-11) 


(6-12) 


If  (12  and  (4)  are  substituted  into  (11)  it  gives; 

P  Cae  -  2^(0  ak)  +  Kc  0e  -  2^(0  uj  ?k) 

-  (O  Jny  -y, 

F  <«m  -  S<pa(0  ac)  +  Kc  0m  -  2<oa(C)  PJ 


=  E<PaCO  <£vy»  <6-13> 

If  C  and  are  specified  by  some  design  requirements,  all  coeffi¬ 
cients  are  fixed  in  (6-13)  and  then  equations  become  simultaneous  equations 

of  P  and  K  ; 

c 


L.  P  +  L,  K  =  N 

X  Z  C  G 


M.  P  +  M„  K 

1  2  C 


K 

m 


(6-14) 


whose  solution  is 


e  2  m  2 

LI  «2  -  L2  «1 

_  Nm  L1  -  Me  M1 
c  “  M2  -  L2 


(6-15) 
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Here 


L1  *  <4v 

L2=  (3e-LWl(O^V 

Mj  =  (  on  -  S<Pa(C),  «£  C<c)  (6-16). 

«2  *  <  Sm  •  =  «>«>  “n  0c> 

Ne=  (  (C)  “4  yk"  ye> 

»„=  (  £<Pa(°  “nyc  *yn) 

If  the  values  found  for  P  and  are  physically  realizable,  then 
the  design  is  completed.  If  it  is  not  the  case,  then  C  or  co  may  be  varied 
to  find  the  acceptable  values,  in  case  such  change  of  values  of  £  or  and 
oun  are  permitted.  If  it  is  not  permitted, “then  multisection  compensation 
is  needed. 

Example  6-1 

420 

Gu(s)  =  s(s  +  l)(s  +  15) 

is  to  be  compensated  to  have  C  >  0.6  by  a  cascade  filter. 

Solution: 


K  (s  +  z) 

Gs 


GC<S>  (s  +  P) 
Characteristic  equation  is; 


s4  +  (16  +  P)  s3  +  (15  +  16  P)  s2  +  (420  K  +  15  P) 

•  c 


s  +  420  P  =  0 
Bj  -  B2  eq uations  for  C  ~  0-6  are: 


B.  =  7“  (420  P)  +  (16  +  P)  CO2  -  1.2  w3 
1  Cu  o 

n 


n 


B„  =  +  1.2  (16  +  P)  CO  -  0.44  co2 

^  (jjL  n  n 

n 

A  =  15  P  +  420  K 
1  c 
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705.7  P  =  7.58 


A2  =  15  +  16  P 


co  =  0.7  is  tried; 
n 


A1  =  Bi  “* 
a2=  B2 

P  =  0.0034 

K  =  0.024 
c 

P  _ 


420K 


+  825  P  =  2J 


0.141 


Thus  the  required  Gc(s)  has  the  form 

G  (S)  0.'024(s  +  0.141) 

(s  +  0.0034) 

The  root  locus  is  seen  in  Fig.  6-2.  The  points  encircled  by  triangles  are 
the  roots  of  the  system.  This  shows  the  applicability  of  the  procedure. 

6.2  Cascade  Compensation  (Continued) 

For  many  problems  it  is  more  convenient  to  use  less  sophisticated 
design  techniques,  such  as  a  combination  of  sketches  with  basic  graphical 
interpretations,  or  curves  calculated  with  the  digital  computer.  Illus¬ 
trations  of  these  techniques  are  given  in  tke  following  section. 

Example  6.2 

G  (s)  = 


106 


Solution; 


u'  '  s  (s  +  5) 

Design  a  lead  compensation  which  makes  the  system  have  £  >  0.6. 

106.Kc(s  +  z) 

G  (S)  =  s(s  +  5)(s  +P) 

The  characteristic  equation  is 


s3  +  (5  +  P)s2  +  (5P  +  106  K  )s  +  106  K  z  =  0. 

c  c 

Mitrovic  equations  for  Bq  -  Bj  scheme  ares  for  this  system; 

Bq  =  (5  +  P)  co2  -  1.2  w3 
B1  =  1.2(5  +  P)  co  -  0.44  co2  . 


Here  (0^(0  are  expressed  numerically  for  the  value  £  =  0»6.  Previous  tables 
tell  that  the  figure  of  Bq  -  B^  is  like  Fig.  6-2  to  have  an  enclosure.  It 
is  actually  seen  that 
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I 


CO  .  <  co 
ml  m2 


holds  for  any  P  by  the  table.  Also  they  show  that 
B. 


1  3 

=  T7T  (5  +  P)  . 
o  max  10  '  7 


If  A  =  106  K  z  =  106  P  is  in  or  on  this  enclosure,  the,  following 
o  c 

relation  must  exist: 

1°6  P  <  Yo“  <5  +  F)3. 

P  >  25,  and  0<  P<  0.1  satisfy  this  inequality.  As  a  load  section  is  re¬ 
quired,  P  =  28  is  adopted.  Then  Aq  =  2968.  Bc  is  equal  to  this  value  if 
co  fa  23.  Substituting  this  value  and  P  =  28  in  the  B^  equation,  it  is  seen 
Bx  =  670 

Aj^  must  have  this  value  if  M  (Aq,  A^)  is  on  the  £  =  0.6  curve. 


which  gives 


and 


A.  =  5  P  +  106  K  =  670 
1  c 


KC=  5* 


z  =  P/K  =  5.6. 
c 

Timas  Gc(s)  turns  out  to  be 

G  (£)  =  5  <»  +  S.6) 
ac' '  (s  +  28) 

If  a  lag  network  is  desired,  P  -  0.1  is  a  good  choice. 

A  =  106  P  =  10.6 
o 

Bq  has  this  value  when  co«  2.  Bj  has  the  value  10.48  correspondingly. 

A.  =  5  P  +  106  K  =  10.48 
l  c 

K  »  0.1 
c 

z  «  P/K  =  1 

Thus  the  lag  compensation  has 


G  (si  =  0-1  (s  4-  1)  « 

6cW  (s  +  0.1) 


Example  6.3 


G  (s)  =  8  (s  +  6) 

u  s(s  +  l)(s  +  3) 

Stabilize  this  system  with  a  single  section  compensation  circuit. 
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Solution 


_  f  v  K  (s  +  z) 

G  (s)  =  c% 

(s  +  p) 

The  characteristic  equation  for  the  whole  system  is 

s4  +  (4  +  P)  s3  +  (3  +  4P  +  8K  )  s2  +  {3P  +  (6  +  z)  SKc  } 

s  +  48  K  z  =  0» 

c 

The  plane  for  £  =  0  is  chosen  and  the  Mitrovie  equations 

are 

Bl  =  (4  +  P)  co2 
B2  =  48  P/co2  +  co2 

K  z  =  P  is  the  relation  to  Maintain  the  error  coefficient  as  before.  From 
c 

the  tables  given  previously,  it  is  seen  the  graph  of  Bj  -  B2  has  the  shape 
shown  in  Fig.  6-3.  Also  the  tables  give 

B1  min=  4««*P)  JTt 

B2  min  =  8‘^r_ 

Polar  analysis,  for  which  explanation  is  already  given,  tells  that  the 
M  (A^,  A2)  point  with  the  following  relation 


lies  in  the  stable  domain. 

These  equations  are  transformed, 

11  P  +  48  K  =  4  (4  +  P)  JW 

c 

3+4P+8K  >8  VsF 

Multiplying  the  inequality  by  6,  subtract  the  former  equation  from  it,  and 
the  following  inequality  is  derived: 

(18  +  3P)  >  748P  (8  -  P) 

P  >  4  and  P  <  0.1  satisfy  this  relation.  Here  it  must  be  noticed  that  too 
large  and  too  small  values  of  P  tend  to  make  the  corresponding  value  of 
Kc  impractical.  In  this  case,  P  =  8  and  P  -  0.1  are  found  to  be  good  trial 
values. 
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> 


(i)  P  =  8 

Blm=  12  X  7384  =  235.2 

Ax  =  88  +  48  Kc  =  235.2 

K  =  =  3.06 

c  48 

A0  =  35  +  8K  =  59.5 

z  c 

B2m=  2  748P  =  39.2  <  A2 

Thus  M  (A^  A2)  surely  lies  its  the  stable  domain.  Gc 
is  £ound  to  be 

r  -  3.06  (s  -t-  2.56) 

Gc(s)  "  (3+8) 

which  is  a  lead  network. 

(ii)  P  =  0.1 

By  similar  calculation^  a  lag  network  with  the  following  trans¬ 
fer  function  is  obtained; 


G  fsx  =  0-17  (s  +  0.61) 

V  '  (s  +  0.1) 


Example  6.4 


G  (s\  _  +  6) _ 

euw  s(s  +  l)(s  +  3) 


is  to  be  compensated  to  have  £  >  0.4  by  a  single  lead  network. 

Solution;  ^ 

B.  =  0.8  +  A,  co2  -  0.8  co3 

1  CO  <5 


B  =  -|  +  0.8  A-co  +  0.36  co2 
2  co 

First,  assume  P  =  15,  as  was  done  in  the  previous  problem.  Comparing 
coefficients 

B.  =  +  19  co2  -  0.8  co3 

1  CO 

720  2 

B  =  —  +  15.2  CO  +  0.36  c 0  . 

2  c C 

These  equations  tell  that  the  B^  -  B2  curve  starts  at  infinity  in  the 
first  quadrant  and  ends  at  infinity  in  the  second  quadrant,  as  is  seen 
in  Fig.  6-4.  From  the  polar  plot  of  Fig.  6-5,  it  is  easily  seen  that 
the  enclosure  exists  only  when  the  B^  -  B2  curve  takes  the  shape  shown 
in  Fig.  6-4.  This  means  B^  is  to  have  one  minimum  and  one  maximum,  while 
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while  has  just  one  minimum.  If  tables  in  THALER  &  RCWNs  Feedback 

Control  Systems  are  used,  this  calculation  is  easily  done  and  the  shape  of 

the  curve  is  shown  as  follows  in  Fig.  6-6.  The  digital  computer  does  the 

same  thing  if  the  attached  programing  is  used  properly. 

Equation  (11)  shows; 

A  =  11  P  +  P  +  48  K 
1  c 

A,  -  3  +  4P  +  8  K  . 

4  C 

This  A^  is  equated  to  =  600  to  give  the  value  of  Kc» 

K  «  9.1 
c 

This  gives  the  value  of  A 2 

A2  =  63  +  72.8  =  135 


which  surely  lies  in  the  £  =  0.4  enclosure.  Thus  the  required  circuit  is 


Example  6-5. 


-  9.1  (s  +  1.65 
c  (s  +  15) 


,  v  420 

GuW  "  s(s -+  i)(s  +  15) 


is  to  be  stabilized  by  a  single  section  compensation  circuit. 
Solution: 

G(s)  =  _ 


420  K  (s.+  z) 
c 


s(s  +  19  (s  +  15)  (s  +  P) 

is  assumed  to  be  the  transfer  function  of  the  stabilized  system.  Then  the 
characteristic  equation  for  the  system  is 

s4  +  (16  +  P)s3  +  (15  +  16  P)s2  +  (15  P  +  420  K  )s 

c 

t 

+  420  K  z  =  0 
c 


As  Aq  is  determined  by  P  only,  the  B,  -  B2  scheme  is  preferable. 
The  tables  give  the  following  equations  for  £  =  0; 


Bj  8  W 

B0  =  A  /cj2  +  CO 2 

2.  O 

and  the  shape  of  the  curve  is  shown  in  Fig.  6-7.  Also  the  tables  tell  that 
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b2«W  * 2  AT 

Bl<“mln)  '  A3  '/ST^ 

The  polar  analysis  indicates  that  the  M  point  with 
A1  =  Bl^°min^  =  B1  min. 


A2  -  B2  min 

certainly  lies  in  the  stable  domain. 

These  equations  can  be  transformed  into 


15  P  +  420  K  =  7420  P  (16  +  P) 

c 

15  +  16  P  >  2  7420  P 
The  inequality  is  satisfied  when 
P  >  4.5  or  P  <  0.2 

From  the  root  locus  considerations  a  large  value  of  P  may  be 
preferred,  but  too  large  a  value  tends  to  make  the  filter  design  im¬ 
practical.  This  is  also  true  with  the  case  of  two  small  a  value  of  P. 
P  =  25  is  quite  a  good  guess  as  is  shown. 


P  =  25 


Bx  =  (16  +  25)  7420  X  25  =  4182 

A  =  15  X  25  +  420  K  =  375  +  420  K 
1  c  c 

K  =  9 
c 

and 

B2  =  2  7420  X  25  =  204 
A2  =  15  +  16  X  25  =  415. 

This  shows  the  M  (A^,  A2)  point  surely  lies  in  the  stable  domain. 

The  G  (s)  required  can  thus  have  the  following  forms 
c 


Example  6.6 


G  (s)  —  9  (s  +  2.77)  . 

(s  +  25) 


Gu<S>  = 


3780  (s  4-  2.77) 

8(8  +  1)(8  +  15) (s  +  25) 


is  to  be  compensated  to  have  £  >  0.6. 
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Solutions 


If  the  compensator  is  assumed  to  be 
K  (s  +  z) 

Gc^s-rrr  > 

the  characteristic  equation  is 

s5  +  (41  +  P)s4  +  (415  +  4lP)s3  +  (375  +  415P  +  3780K  )s2 

c 

+  (4155P  +  10500K  )s  +  10500P  =  0 

C 

The  B^  -  B2  equations  are  (for  C  =  0.6) 

B  =  1.2  A  /t0  +  A,co2  -  1.2  A. co3  +  0.44  to4 

1  O  3  4 

B«  =  A  /to2  +1.2  AM  -  0.44  AM2  -  0.672  td3. 

Z  O  J  4 

These  equations  are  too  complicated  and  it  is  rather  unwise  to  have 
preliminary  checks  such  as  to  see  what  ranges  of  P  are  adequate,  or 
whether  one  section  will  suffice.  It  is  noticed  that  a  larger  value 
of  P  is  needed  for  the  purpose  from  the  view- point  of  root  locus  method. 
P  =  30  is  tried  and  the  graph  is  depicted  as  in  Fig.  6-8. 

to  2.5  5  7.5  10  12.5  15 

B1  16  x  104  11  x  104  11  x  104  12  x  104  13  x  104  12.8  x  104 

B2  5.4  x  104  2.2  x  104  2  x  104  1.9  x  104  2  x  104  2.1  x  104 

This  P  =  30  is  a  little  short  of  the  specification  as  the  K£  line  can'°t 
hit  the  enclosure.  The  value  of  P  must  be  changed  a  little.  So  the 
effect  of  changing  P  is  ckecked: 


d  B1  _ 

1.2  x  10500 

+  41  ,to2 

-  1.2  td3 

dP 

CO 

d  b2  _ 

10500  +  U2 

x  41  to  - 

0.44  to2 

dP 

J 

This  shows  that  at  to  =  10 

Bj  (P  +  A  P)  -  Bl(P)  =  15500  A  P 
B2  (P  +  A  P)  -  B2(P)  =  157  A  P 
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Also  it  is  seen  that 


A.  (P  +  A  P,  K  +  A  K  )  -  A,  (P,  K.  )  =  4155  A  P  +  10500  A 

1  C  C  1  c  c 

A2(P  +  A  P,  Kc  +  A  Kc)  -  A2  (P,  Kc)  =  415  A  P  +  3780  A  K£ 

If  K  =  lo63,  then 
c 

A1(P,  Kc)  =  Al(30,  1.63)  =  141765 
A2(P,  Kc)  =  A2(30,  1.63)  =  19000 
Thus  it  is  seen: 

B^P  +  A  P)  =  12  x  104  +  4160  A  P 

A.(P+AP,K  +  A  K  )  =  141765  +  4155  A  P  +  10500  A  K 

X  C  C  C 

A2(P  +  A  P)  =  1.9  x  104  +  157  A  P 

A„(P  +  AP,  K  +  A  K  )  =  19000  +  415  A  P  +  3780  A  K 
z  c  c  c 

Here  (P  +  A  P)  and  b2  (P  +  A  P)  mean  E,  and  B2  on  the  new  curve. 

To  have  the  M  (A^,  A2)  point  in  the  enclosure  the  following 

equations  are  needed s 

A.  (P  +  A  P}  K  +  A  K  )  =  B,  (P  +  A  P) 

X  C  C  X 

A2  (P  +  A  P,  Kc  +  A  Kc)  =  B2(P  +  A  P) 

Substituting,,  the  f  ollowing  results  are  obtained: 

21765  -  11345  A  P  +  10500  A  K  =  0 

c 

258  A  P  -  3780  A  K  =  0. 

c 


If  AK£  =  -0.1,  then 
A  P  =  2 


is  obtained  from  the  first  of  these  equations. 

side  of  the  second  equation  equal  to  894.  This 

lies  inside  the  new  B^  -  B2  curve  by  the  amount 

stability  with  £  —  0.6,  as  the  numerical  values 

Thus,  it  is  seen  that 

=  1.53  (s  +  20,9) 

VS)  (s  +  32) 

is  the  required  compensation. 


This  makes  the  left  kand 
means  the  MCA, ,  A2>  point 
of  894,  which  assures  the 
of  the  graph  shows. 
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Note:  This  method  is  applicable  only  when  A  P  doesn't  change 
the  shape  of  the  resulting  graph  too  much,  that  is£  the  relative  positions 
of  maxima  and  minima  must  be  kept  on  before  and  after  the  change  of  P. 

6.3  Feedback  Compensation 

Mitrovic's  method,  whether  it  implies  the  original  or  extended 
one,  is  primarily  just  concerned  with  the  system  characteristic  equations 
and  so  general  procedures  are  the  same  in  both  cascade  and  feedback  com¬ 
pensation.  But  in  the  latter  compensation  if  a  pure  zero  filter  or  pure 
derivative,  and  accelerative  ones  are  permitted,  the  order  of  the  orig¬ 
inal  system  need  not  be  heightened.  This  means  seme  simplification. 

Example  6.7 

For  the  system  of  Fig.  6-9,  find  an  adequate  value  for  to  make 

the  system  have  £  >  0.6. 

Solution:  The  system  characteristic  equation  is 

s3  +  16  s2  +  (15  +  420  K  )s  +  420  =  0 
Bq  -  B^  equations  with  £  =  0.6  are  given  to  be 

B  =  16  co2  -  1.2  u ? 
o 

B^  =  19.2  co  -  0.44  co2 

The  polar  analysis  shows  that  B  -  E,  curve  must  have  the  form  of 
Fig.  6-10  to  satisfy  the  specification.  The  values  of  co^  are  given  by  the 
tables:  1 

tomi  =  8.8 

^m2  =  22 

which  show  the  relation  among  the  three  coefficients  is  favorable  to  have 
an  enclosure. 

B  (co  .)  =  668.3. 
o  mi 

As  Ac  is  kept  to  be  420,  there  must  be  two  points  on  the  curve  for  this  case, 
co  =  5.9  is  found  to  satisfy  Bq  =  420.  B^(co)  for  this  value  is 

Bx(5.9)  =  19.2  x  5.9  -  0.44  x  (5-9) 2 
=  98 

Ax  =  15  +  420  Kt  =  98 
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Thus  Kt  is  found  to  be; 

Kt  «  0.2 

Example  6.8 


a  M  * _ 4SQ.feJi_» _ 

Vs'1  s(s  +  l)(s  +  5)(s  +  7) 

is  stabilized  by  acceleration  feedback,  as  shown  in  Fig.  6-11 
What  values  of  are  applicable? 

Solution: 


The  system  characteristic  equation  is 

4  +  (13  +  45 OK  )s3  +  (47  +  1800  K  )s2  +  485  s  +  1800  =  0 
s  a  3 

“  Bj  plane  is  preferable  in  this  example,  which  gives  for  £  =  0, 

B2  =  -j  +  co2 
£0 

b3=  tl 
J  2 
£0 

From  the  table  it  is  seen  that: 

2 


and 


“Sln='/l800i=  42.-4" 


=  84‘8 
83<“»in)  *  11 -4 


From  the  polar  analysis  of  stability,  it  is  seen  that  the 
curve  must  have  the  shape  shown  in  Fig.  6-12,  and  M  (A^»  A^)  must  be  inside 
this  curve.  So  the  following  calculations  are  made  to  find  the  value  of  K 


a 


2 

£0 

1 

810 

42.4 

50 

100 

200 

18Q0 

2 

CO 

1800 

180 

42.4 

36 

18 

9 

B2 

1801 

190 

84.8 

86 

118 

209 

B3 

485 

48.5 

11.4 

9.7 

4.85 

2.43 

The  graph  of  Fig.  6-13  shows 


K  > 


17 


=  0.04 


a  —  450 

is  applicable  to  get  a  stabilized  system. 
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U)=oo 


60=0 


A 


Figure  6-4 


Figure  6-5 


Figure  6-6 


o  o 


PROGRAM  MIT  BLM 

PROGRAM  IS  TO  GRAPH  MITROVIC  THALER  BL  VS.  BM  CURVES  FOR  EQUATION 
A(N)*S**ty-A(N-l)*  S**(N-1)+  . .  +  BL*S**’L  +  .  o+BM*S**Mf  .  .+A(1)* 

S  +  AZERO3 C, 

C  ACCORDING  TO  METHOD  DEVELOPED  BY  THALER  AND  OTHERS. 

C  PROGRAM  IS  DIMENSIONED  FOR  MAXIMUM  N  =  30. 

C  PROGRAM  IS  NORMALIZED  SO  AS  A(N)  =1.0 

C  A(L)  AND  A(M)  ARE  ASSUMED  TO  EXIST,  BOTH  BEING  EQUAL  TO  ZERO. 

C  DATA  CARDS 

C  1. USERS  NAME  AND  PROGRAM  IDENTIFICATION  COLUMNS  1  THRU  27. 

C  2. OUTPUT  SPECIFICATION 

C  BLANK  CARD  IS  PRINT  AND  GRAPH  ARE  BOTH  REQUIRED. 

C  WORD  GRAPH  IN  COLUMNS  1  THRU  5  IF  GRAPH  ONLY  IS  REQUIRED. 

C  WORD  PRINT  IN  COLUMNS  1  THRU  5  IF  PRINT  ONLY  IS  REQUIRED. 

C  3. ORDER  OF  EQUATION,  N,  12  FORMAT. 

C  4.  COEFFICIENTS  AZERO  THRU  A(N),  8F10.2  FORMAT 

C  5. VALUES  OF  L  AND  M,  212  FORMAT. 

DIMENSION  A(30) ,  PHI(60),  BLAPHI  (60),  BMAPHI(60),  BL(900) ,  BM(900), 
1  ITITLE(IO),  JTITLE(IO);  B(31)  > 

C  READ  INPUT  DATA 

READ  100, (ITITLE(J) ,  J=l,  4) 

READ  100, 1  OUT  PUT 
READ  101, N 

READ  102,  AZERO,  (A(J),J  =  1,N) 

READ  103, L,M 

100  FORMAT  (A7,3A8) 

101  FORMAT(12) 

102  FORMAT (8F 10. 2) 

103  FORMAT(212) 

C  SET  UP  ITITLE  FOR  GRAPH  AND  JTITLE  FOR  PRINT. 

JTITLE(l)  ,=8HM 

LDQ (ITITLE  +  1),ENA(20b),  LRS(6),  STQ (ITITLE  +  1) 

LDA (ITITLE  +  4),  ARS (6) , LDQ (JTITLE  +  1) ,LLS(6) ,STA(ITITLE  +  4) 

ITITLE (5)  =  8HITROVIC 

ITITLE (6)  =8HBL  VS  B 

ITITLE (7)  =8HM  CURVE 

DO  2  J=8,  10 

2  ITITLE (J)  =8k 
DO  3  J  =1,10 

3  JTITLE (J)  =  ITITLE (J) 

C  WE  CHECK  OUTPUT  REQUIREMENTS 

J-TESTl  =  8H 
ITEST2  =  8HGRAPH 
TTEST3  =  8 H PRINT 
IF(ITESTl IOUTPUT)12, 11,12 

11  INDICPR  =  1 
INDICGR  =  1 
GO  TO  17 

12  IF(ITEST2  -  IOUTPUT)14,13,14 

13  INDICPR  =0 
INDICGR  =1 
GO  TO  17 


6-18 


o  o 


14  TF(ITEST3-I0UTPUT) 16,15,16 

15  TNDTCPR  =  1 

TNDICGR  =  0 
GO  TO  17 

16  PRINT  208 

STOP 

208  FORMAT  (20H  ERROR  IN  DATA  CARDS) 

C  PRINT  RECORD  OF  INPUT  DATA. 

17  PRINT  200 

PRINT  204,(JTITLE<fJ),J  =  1,7) 

PRINT  206 

PRINT  207 ,AZERO, (A(J) , J  =1,N) 

PRINT  20 2, L 
PRINT  203, M 

200  FORMAT (lHl) 

201  FORMAT (IHO) 

204  FORMAT (4X.7A8) 

205  FORMAT (1H0S25HTHE  ORDER  OF  EQUATION  IS  ,12) 

206  F0RMAT(1H0,39HTHE  COEFFICIENTS sAZERO  THRU  A(N)  ARE  ) 

207  FORMAT (20X.E14. 6) 

202  FORMAT (7HBL  IS  Bs12) 

203  FORMAT (7HBM  IS  B,12) 

C  THE  EQUATION  IS  NORMALIZED  SO  AS  TO  GET  ACN)  =  1.0 

ANRECTP  =  1.0/A(N) 

PRINT  212 

212  FORMAT (flHO,  10X.45  HTHE  RESULTS  REFER  TO  THE  NORMALIZED  EQUATION. ) 

AZERO  =  AZERO*ANRECIF 
DO  5  J=1,N 

5  AfJ)  =  A(J)*ANRECIP 

THE  STEP  SIZE  WILL  BE  TAKEN  AS  O.OOl  TIMES  THE  FREQUENCY  OF  INTER¬ 
EST.  FOR  THIS  AIM  WE  LOOK  FOR  THE  VALUES  OF  OMEGAN  WHICH  MAKE  BOTH 
OF  THE  LAST  TERMS  TN  BL  AND  BM  DOMINATE  FOR  EACH  VALUE  OF  ZETA. 

LM  =L-M 
LMl  =XABSF(CLM) 

IF((LMl/2)*2-LMl)175,74,75 

74  ZETA  =0.1 
KK  =5 
GO  TO  76 

75  ZETA  =0.0 
KK  =6 

76  DO  50  JJ  =1 ,KK 

C  CALCULATE  CFBYSHEV  POLYNOMIALS  e 

20  PHK31)  =-l. 

PKI<f32)  =2.*ZETA 
NP30  =  Nf 30 

DO  21  K  =33,  NP3C 

21  PHI (K)  = - 2 . 0*ZETA*PHI f K- 1) - PHI (K- 2 
PHI (30)  =0. 

UN  =30-N 
DO  105  J  =MN,29 
105  PH If J )  =-PHI(60-J) 
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C  THE  TERMS  IN  BL  AND  BM  ARE  CALCULATED  IN  TWO  STAGES  „ 

DO  104  1  =1,N 
104  Bfl  +  1)  =  A(l) 

B(l)  =  AZERO 
NPI  =  N  +1 
NMl  =  N-  1 
NM2  =  N-2 
L29  =  29-L 
M29  *  29-M 
DO  22  1  =1,  NP1 

BLAPHia)  »  BCI)  *PHIfI/+  M29 ) / PHI ( 30- LM) 

22  SMAFHI(X)  =  B(I)  *PHI(I  +  L29)/PHl(30fLM) 

OMEGAM  =0.5 

DO  8K  =  1,12 

OMEGAM  =  2.0*0MEGAM 

SL  =  BLAPHI(l)  *  OMEGAM** (  -L) 

Ll  =  L  +  1 
DO  311  1  =  2.NP1 

311  SL  =  SL  +  BLAPHI(I)*OMEGAM**(T  -Ll) 

NL  =N-L 

IF CBLAPHI(NPl) >71,170,71 

170  IF(BLAPHI(N) >70,172,70 

71.  BSA  =  BLAFHI (NPI )  *OMEGAM**NL 
GO  TO  72 

70  BSA  =  BLAPHI(N)*OMEGAM**(NL  -  1) 

GO  TO  72 

172  mBLAPHI(NMl))171,173,171 

171  BSA  =  BLAPHI (NMl)*OMEGAM**(NL  -2> 

GO  TO  72 

173  BSA  =  BIAPHI(NM2)*OMEGAM**(NL  -  3) 

72  BSLN  =  ABSF(BSA) 

BSL  = ABSF ( SL- BSA) 

TF(BSL-BSLN)9,8,B 

8  CONTINUE 

9  STEP  =  OMEGAM 
GMBGAN  =0.5 
DO  68  K  =1,12 
OMEGAN  =  OMEGAN*  2 . 0 

SM  =BMAPHI(1)*0MEGAN**(-M> 

Ml  =  M  +  1 
DO  312  1  =2,NPl 

312  SM  =  AM  +BMAPHI(fI)  *  OMEGAN**  (I -Ml) 

IF ( BMAPHI ( NPI ) ) 81 , 180 , 81 

180  IF  (BMAPHI  (N»  80,  182,80 

81  BSB  =  BMAPHI (NPI )*OMEGAN**  NM 
GO  TO  82 

80  BSB  =  BMAPHI (N)*OMEGAN**(NM-l) 

GO  TO  82 

182  IF (BMAPHI (NMl) >181,183, 181 

181  BSB  =  BMAPHI (NMl )*OMEGAN**(NM- 2) 

GO  TO  82 

183  BSB  = BMAPHI (NM2)*0MEGAN**(NM-3) 
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82  BSMN  =  AE3F  (  BSB) 

BSM  =  ABSF  (SM  -  BSB) 

IF  (BSM  -  ESMN) 69,68,68 

68  CONTINUE 

69  STEM  =CMEGAN 

IF (STEP  SIEN)216,215 ,215 

2X6  STEP  =  STEN*0.  OOl 
GO  TO  123 

215  STEP  =STEF*O.OQl 
123  DO  310  IOMEGAN  =1,900 
OMEGAN  = IOMEGAN 
(MEGAN  = OMEGAN*  STEP 
SUML  =  BLAPH1  (1)  *OMEGAN**(-L) 

SUMM  =BMAFHI  (1)  ♦OMEGAN*'*  (  -M) 

DO  23  1=2, NP1 

SUML?  SUML  +BLAPHI  Cl)*  0MEGAN**(I-L1) 

23  SIW3UMM  +  BMAPEI  < I  )*OMEGAN**  ( I  -  Ml ) 

BL(  IOMEGAN)  =  SUML 

BM( IOMEGAN)  =  SUMM 

C  PRINT  EVERY  TENTH  POINT  IF  PRINT  OUT  IS  REQUIRED. 

IF(INDICPR-I)  310,24,310 

24  IF ((MEGAN-1)  25,27,25 

25  IF (XMODF (LOMEGAN  ,  100»26,29,26 

26  IF (XMODF (IOMEGAN  ,  10))310,30,310 

27  PRINT  200 
PRINT  209 , ZETA 
PRINT  210 

GO  TO  310 

29  PRINT  201 

30  PRINT  211,  (MEGAN , BL(IOMEGAN) , BM (I OMEGAN) 

209  FORMAT (10X, 35HRESULTS  OF  COMPUTATION  WITH  ZETA  =  ,F4.2) 

210  FORMAT (//20X,  „10H  OMEGAN  ,10H  BL  ,10H  BM  //) 

211  FORMAT  (19X,  Ell.  5 ,2X,Ell  »5,2X;,E11 .5) 

CONTINUE 

NUMPTS  =  900 

C  GRAPE  IF  GRAPH  IS  REQUIRED 

IFQMDICGR  -1)49,132,49 
132  IF{ (LMl/ 2 )*2-LMl ) 134 , 232 , 134 

134  GO  TO  (33,34,35,36,37,38) ,JJ 

232  GO  TO  (33, 34, 35, 36,37), JJ 

33  MODCURV  =0 
LABEL  =  2H1 
GO  TO  39 

34  MODCURV  =0 
LABEL  =2H2 
GO  TO  39 

35  MODCURV  =0 
LABEL  =  2H3 
GO  TO  39 

36  MODCURV  =0 
LABEL  =  2H4 
GO  TO  39 

37  MODCURV  =0 
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LABEL  -2H5 
GO  TO  39 

38  MODCURV  =0 
LABEL  =2h6 

39  SFX  =0.0 
SFY  =0.0 
MINOFFX  =0 
MINOFFY  =0 
LABELNO  =11 
MODE  =  0 

Nl  =0 
ME  =0 

CALL  GRAFH2  (NUMPTS , BL , BM,  8 , MODCURV , LABEL , ITITLE , SFX , SFY , MINOFFX , 
1MIM0FF Y , LABELNO , MODE s  Nl , N2) 

49  ZETA  =ZETA  +0.2 

50  CONTINUE 
END 

END 
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APPENDIX 


If  cor  goes  to  zero  in  the  equation  (3-5), 


it  becomes 


?H>  “  <Ai  ‘  V  "n  *t  +  iAJ  •  Bj>  “4  ‘j 

ies 

lira  F  (w  )  =  -lin  { B.  (0/  e.  +  B.(w  )  J  e) 
w  _0  n  li'n  1  T  j  n  n  jJ 


n 


.  =  A 


<P  AO  2.  -  ^iiO  ^ 


3-jP  1  V)(°  i 


■  A°  U 

,  ”i-j(D  L 


j(D  et  -<»,«>  ej 


Here  e.  and  e.  are 
i  J 


et  =/  i  f  +  ie  -  ( -  c  +  i  fi  <2  )* 

e j  =  /j  \  +  je  =  (  -  C  +  i  J 1  -cV 

and  the  following  expressions  are  derived 


(3-5) 


(A-1) 


(A-2) 


=  <  + 1  y  1  -  cz »  <p0«)  =  c  <px(o  - 1  yr^^co 

e2  =  -1  +  2C2  +  I  VT-  ^ (O  =  *!«>  +  C  <P2(C)  “ 

-i  7i  -cz<p2co 

(  1  M  I 

®n  =  Vl  i0  +  C  ^n(C)  "  1  ^  -C2<Pn(C).  (A-3) 

From  these  formula,  it  is  seen  that 


<Pt(C)  e.  =  pt(0  +  C  ^(C)  <Pj(0  - 

-  l7l  -^V^OVA-O 

^(C)  e.  =  ^(O  <Pi.1(C)  +  C  ^(C)  <Pi(C)  - 

-  1/1  -c2  <Pj<C)  ^(D 


(A-4) 


APPENDIX  (Continued) 

Subtracting  the  latter  from  the  former,  it  is  seen 

e.  -  <Pj(C)  e1=  <pi(C)  -  <Pj(C)  <Piml(0  (A- 7) 

Also  it  is  seen  that 

*i<C>  Oj.jCC)  -  ^(C)  *^<0  *  PtmjiO  (a- 8) 

from  a  simple  mathematical  induction. 

Substituting  (A-7)  and  (A-8)  into  (A-l)  gives 

lim  F  («  )  =  A 

ui  -*  0  ”  ° 

n 

which  is  the  required  result,  equivalent  to  equation  (3-6) 


TOHOMI  OHTA  was  born  in  Salisbury  Avenue,  Kaoloon,  Hongkong 
in  1923.  He  was  educated  at  Keio  University,  Tokyo,  receiving  his 
bachelor's  degree  from  the  Electrical  Department,  Faculty  of  Engineer¬ 
ing  in  1948,  and  later  a  Master's  Degree.  He  has  been  an  Instructor 
in  Physics  at  Keio  University  (1953  -  1957),  and  joined  the  Electrical 
Department  of  the  Defense  Academy  of  Japan  as  an  Assistant  Professor 
in" 1957.  The  research  reported  here  was  done  while  visiting  The  U.S. 
Naval  Postgraduate  School  in  1963.  Professor  Ohta  has  also  published 
papers  on  "Electric  Discharge  in  Gases"  (M.  S.  Thesis,  Keio  University), 
'Nuclear  Shell  Structure"  (Physical  Society  of  Japan,  and  Memoirs  of 
the  Faculty  of  Engineering,  Keiv  University),  "On  the  Design  of  a 
Homogeneous  Nuclear  Reactor"  and  "On  the  Design  of  a  Heterogeneous 
Nuclear  Reactor"  (Memories  of  the  Defense  Academy  of  Japan). 


GEORGE  J.  THALER  was  born  in  Baltimore,  Maryland  in  1918.  He 
was  educated  at  The  Johns  Hopkins  University,  receiving  degrees  of 
Bachelor  of  Engineering  (1940)  and  Doctor  of  Engineering  (1947).  He 
was  an  Instructor  at  Johns  Hopkins  (1942-1947),  an  Assistant  Professor 
at  the  University  of  Notre  Dame  (1947-1951)  and  has  been  on  the  Faculty 
of  the  U.S.  Naval  Postgraduate  School  since  1951.  He  has  specialized 
in  Automatic  Feedback  Control  Theory,  and  has  written  four  textbooks 
and  over  thirty  research  papers  in  this  area. 


% 


